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Section 2.6 


All temperature variances quoted in this section should read 
(°c? x 103). 


Section 3.3 


Pressures plotted in figures 25a - 25q are in millibars (mb). 
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PREFACE 

The object of this report is to communicate certain aspects of the 
oceanographic observations taken along, and in the vicinity of, Line P during 
May 1972. For this reason, some effort has gone into reducing and inter- 
preting the data obtained. It is hopal that the results presented will 
facilitate a better understanding of the physical oceanography in the 
oceanic region adjacent to the coast of British Columbia and also aid 
in the planning of similar surveys in the future. 

As with any work of this type, there has been a large amount of 
input from other individuals. For ereenetih collecting the data, I 
would like to thank the officers and crews of the weatherships CCGS 
VANCOUVER and QUADRA with special thanks to Captain A.A.R. Dykes of the 
QUADRA for his cooperation in allowing a replacement sensor to be trans- 
ferred aboard his ship while at sea and for occupying all inbound Line P 
stations. In this connection also, the effort of Mr. W. Hansen in his 
role as the sea-going technician on the QUADRA is gratefully acknowledged. 

The amiable and generous assistance of Captain M. Dyer, his 
officers and crew of the CFAV LAYMORE made for both a useful and pleasurable 
cruise. My co-worker, Mr. C. de Jong is especially thanked. Without his 
expertise and hard work much less would have been accomplished. 

With regards the onshore reduction of the data, I would like to 
thank Mr. B. Minkley, Ms. M. Dyer and Ms. R. Thomas. Mssrs. D. Smith, 
E. Luscombe, B. Johns and A. Douglas and Ms. L. Kuwahara are especially 
acknowledged for their excellent work in the computational portion of the 


data analysis. Moreover, the assistance of Ms. Kuwahara, Ms. S. Greckol 


and Mr. P. Huggett in plotting the results has been very much appreciated, 
as has that of the photographic section of the Hydrographic Service 
(Victoria). Finally my thanks to Drs. J. Garrett and S. Tabata for 


reviewing the manuscript and for making useful suggestions. 
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1. INTRODUCTION 

Canadian operation of Ocean Weather Station P (latitude 50°00' N, 
longitude 145°00' W) commenced in December 1950 with the naval frigate class 
ships CCGS St. CATHARINES and CCGS STONETOWN alternatively occupying the 
station for periods of six weeks at a time. A program of twice-daily bathy- 
thermograph (BT) observations was added to the ships' routines in July 1952. 
In August 1956, St. CATHARINES further augmented its oceanographic observ- 
ations with a more extensive BT program which by April 1959 had been extended 
to include a series of "Line P" stations en route to and from Station P. 

The present weatherships, CCGS VANCOUVER and CCGS QUADRA, which 
replaced the frigates in April 1967 and October 1967, respectively, now 
occupy twelve positions along the Line (see figure 1). In addition, these 
ships have since augmerrted the BT observations with a program of Expendable 
BT (XBT) and Salinity-Temperature-Depth (STD) observations. The latter 
program, which is still in operation, began in May 1968 onboard the VANCOUVER 
and in February 1969 onboard the QUADRA, using Bissett-Berman model 9006 
or 9040 STD's. 

The results of this report stemmed from a desire to make use of 
the weathership observations along Line P to gain some understanding of the 
“spatial distribution of the physical oceanographic properties off the coast 
of B.C. and their temporal variability over periods of days. As a conse- 
quence of the ships' routines and speeds, however, weathership observational 
data alone is not sufficient to give the required information. The difficulty 
is that of the six weeks spent at sea by each ship, only about two and one- 
half days are required to complete all Line P stations during the outbound 
or inbound leg of a particular cruise. This, together with the requirement 


of keeping P continuously occupied leads to a difference of only a few hours 
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between observations at the most oceanic station (12) but five days between 
observations at the most coastal station (1). To obtain temporal variations 
based on more evenly distributed time lags at all stations, it was therefore 
necessary to employ the services of a third vessel. By arranging the depart- 
ure time of this vessel to P to be a few days in advance of the outbound 
weathership, and its arrival time in port to be nearly coincident with the 
inbound weathership, the desired lags (days) were achieved (see figure 2). 

Although it would have been preferable to have a ship capable of a 
Speed comparable to that of the weathenshine (16 knots), only the relatively 
slow (11 knots) CFAV LAYMORE was available. Nonetheless, the disadvantage 
in speed was somewhat compensated by the LAYMORE's maneuverability and 
working area arrangement which permitted stations to be one rather efficiently. 
The month of May was chosen for two reasons. The first and most important 
consideration was the weather which at that time of year is usually pleasant 
enough to permit all Line P stations to be done by the weatherships; in high 
winds their large amount of superstructure makes observing very difficult. 
Second, May is a time midway between the Enaina Stability in the upper 
Oceanic layers that occurs in early March and the maximal Stability in 
these waters that occurs in late July. The depth of the surface mixed layer is 
for example, shallow enough to be distinguishable from more permanent thermo- 
cline(s) but deep enough not to be affected by the local short period (hours) 
surface heating. 

As with the weatherships, observations from the LAYMORE during the 
May 8-18 cruise were taken using XBT's (Sippican Corporation) and a Bissett- 
Berman STD (model 9006 in this case). The latter measures temperature via a 
thermistor, conductivity through a temperature compensated conductivity cell 


and depth via a frequency modulated transducer Signal. Output is recorded 
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as an analog trace of temperature and Salinity versus depth on lined chart 
paper. These are then digitized ashore using a digitizing table. As the 
accuracy Of the STD's is * 0.03°C for temperature, * 0.05°/00 for salinity and 
t 1% for depth, there should in theory be no loss of accuracy during digitizing 
Since the charts can be read to higher precision than this. In practice, 
however, the fact that the salinity traces 'spike' each time there is a 
Sharp temperature gradient leads to a Significant information loss in the 
case of salinity. This feature of the Bissett-Berman's, arising through the 
difference in the time constants between the temperature sensor and the 
compensated conductivity sensor, necessitates a hand smoothing of the salinity 
traces and therefore to a loss of detail. There is also a tendency for all B.B. 
sensors to drift with time. For each cast, then, a correction factor should 
be determined through a simultaneous Nansen bottle cast. As this is extremely 
time consuming, a few Nansen casts only are usually taken to give an indic- 
ation of the mean correction factors to be applied. On the LAYMORE cruise, 
surface bucket samples were also taken at the time of STD stations. These 
had the advantage that they gave a truer indication of the temperature and 
Salinity of the top layer of water than possible with the STD. The reason 
for this, of course, is that the STD being subject to the rollof the ship 
had to be commenced with the instrument a meter or so below the surface. In 
the data, there has been no attempt to correct for this since it is virtually 
impossible to visually ascertain the actual depth of the STD once it is in the 
water. 

For the purposes of this report, the results of the data reduction 
are presented in two sections. The first deals basically with the spatial 
distribution of physical properties in the Vicinity of Line P although there 


is also some discussion of time variability. In the second section, results 


from all three ships are compared to give some indication of the temporal 
changes that occurred. In both cases, explanations are presented where 


warranted, and possible. Digitized plots of the original data are given 


in the Appendix. 


2. RESULTS from the CFAV LAYMORE cruise: May 8-18, 1972 
2.1. General Discussion 

LAYMORE outbound 

The LAYMORE departed Esquimalt, B.C. at about 1800 GMT on May 8 
bound for Ocean Weather Station P. Station 1 was occupied at about 0200 GMT 
the next day. Except for the continental-margin stations 1, 2 and 3 which 
were taken to 100, 100 and 1000 m, respectively, all STD casts went to 1500 m 
with station 4 further extended to 2400 m using Nansen bottles (at 1200, 
1400, 1600, 1800, 2000 and 2400 m). At each STD station a surface bucket 
Sample was obtained, the temperature read using a thermometer and a sample 
bottled for salinity determination ashore. In addition, Expendable Bathy- 
thermographs (XTB's) were taken at each Line P outbound station as well as 
at intermediate positions (figure 3). An attempt to extend station 9 to 
2400 m using Nansen bottles had to be abandoned because of bad weather. The 
outbound leg of the cruise was completed with a rendez-vous with the QUADRA 
at the edge of the Station P grid at about 0500 GMT on the 12th. A re- 
placement STD sensor for the malfunctioning unit on the weathership was then 
successfully transferred from the LAYMORE and installed. 

LAYMORE inbound 

The inbound portion of the LAYMORE cruise was altered somewhat from 
the outbound track so that some information on the water properties north 
and south of the Line could be obtained (figure 3). North-south legs of this 
course which intersected Line P so as to extend 14° of latitude on either 
side of it, were joined by legs parallel to the Line. The longitudinal 
positioning of stations in this case were determined by the longitude of 
the Line P stations and the latitudinal positions by the requirement of a 


Spacing similar to those along the Line. 
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As with the outbound part of the cruise, a bucket sample was taken 
at each STD station. XBT's, however, were only taken between stations because 
of the limited supply of probes resulting from the fairly high rate of 
failure (~ 20%). Also, an attempt to extend station 18 with a Nansen cast 
had to be abandoned at 300 m due to the large amount of wind drift and 
subsequent large wire angle (~70°). Moreover, station 1 inbound was cancelled 
because of high seas. 

Corrections to STD readings 

Despite the fact that igta for only one shallow and one deep 
Nansen cast are available from the LAYMORE survey, it has been possible to 
determine mean corrections to the STD readings using these and the greater 
number of comparisons performed with the same instrument during the June 23- 
August 10 patrol (72-005) of the QUADRA. The latter are plotted in figures 
4a, b as Nansen readings minus STD readings versus depth. The heavy solid 
line on each curve is the mean correction to be applied for that particular 
cruise. Since the points from the 0-300 m Nansen cast of the LAYMORE fall 
within the range of those from the QUADRA, the same surface error is assumed 
for both cruises. The slight deviation with depth between the corrections 
applied to the LAYMORE observations and those applied to the QUADRA observ- 
ations (figures 5a, b) are a result of assuming that near 1400 m the points 
from the former are more correct than those from the latter. Strictly 
speaking, the difference at this depth is rather insignificant when compared 
to the cast-to-cast variations shown by figure 4. Thus the QUADRA corrections 
could just as readily be used. Moreover, comparison with an even later cruise 
(QUADRA, 15 September-2 November), in which the same instrument was used, 
shows only a slight fertation between mean curves, although again the cast- 


to-cast variations are large. Finally, in accordance with the present 
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procedure followed by the Weathership Group, spikes in the salinity traces 
have been removed by hand before digitizing. A salinity value has been consi- 
dered to be a spike when it is associated with a sharp, localized temperature 
gradient provided the initial change in salinity is in the same senseas the 
initial temperature change and is then followed by a Salinity change in the 


opposite sense. 
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FIGURE 4a: Nansen bottle salinity minus S.T.D. salinity: Solid line (—) 
srepresents mean correction curve for "QUADRA" cruise 72-05 
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Nansen bottle temperature minus S.T.D. temperature: Solid line (—) 
represents mean correction curve for QUADRA cruise 72 - O05: 
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FIGURE 5a: Nansen bottle temperature minus S.T-D. temperature for "LAYMORE" cruise 72-10- 
x surface waters,@ deep waters: Correction curves: 
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FIGURE 5b: Nansen bottle minus STD salinity for"LAYMORE" cruise 72-10; 
x surface water, o deep water; correction curves: 
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2.2. The Spatial Distribution of Physical Properties 

Use of the STD and XBT data makes it possible to obtain the mean 
horizontal distribution of water properties in the vicinity of Line P during 
early May 1972. For consistency, outbound values have been used where there 
are two values for a given station. Plots presented are for a number of 
standard depths over the 0-1500 m range. 

As with any presentation of this kind, only an estimate of the 
mean horizontal distribution of a particular property is obtainable since 
variations due to internal tides, internal waves, inertial waves, horizontal 
advection, etc. cannot be removed. The effects of vertical motions are 
especially strong in regions of large vertical gradients such as occur 
in the upper 100-200 m of the subarctic Pacific. An example of this is 
figure 4 which shows the large differences between STD and Nansen observ- 
ations taken in the highly structured upper layer at slightly different 
times at the same location (Station P). In addition, near-coastal water 
properties can be significantly changed by tidal mixing and river discharge. 
Only below the main pycnocline, (> 200 m) where temporal variations should 
be relatively small, would a single cast be expected to closely reproduce 
the mean values. Even theione is limited by the accuracy of the instrument 
since with increasing depth below the pycnocline observational errors begin 
to become more comparable to changes that can be produced by advection. 
These remarks are partly demonstrated in Table 1 which, as an example ,gives 
the differences in temperature at selected depths between the two different 
samplings of six Line P stations. 

In summary, then, any plot of the horizontal distribution of a 
water property over an extensive oceanic region will contain the integrated 


effect of temporal variations over the observational period. Confidence can, 
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therefore, only be given to the gross features of the observed distributions. 

Temperature Sections (STD) 

The horizontal distribution of temperature at depths obtained from 
the STD observations is presented in figures 6a-i. Except at 1500 m, where 
there is little structure, the isotherms over most of the region show a 
northward decrease in temperature with the nearly zonal orientation in the 
interior tending more north-south as the coast iS approached. The main 
deviation from this trend occurs in the vicinity of consecutive station 26 
(50°47' N, 132940' W) and appears to extend from about 100 to 500 m. 

When the isotherm pattern at Om is compared with that at other 
depths, the time integrated effect of surface heating over the observation 
period becomes apparent. In particular, the change in orientation of the 0m 
isotherms that occurs along Line P between 128°W and 136°W has probably been 
exaggerated as a result of the time delay between Line P outbound observations 
and those to the north and south. At subsurface depths, where the heating 
trend in May has not been felt, the change in isotherm orientation is more 
gradual and undoubtedly gives a more accurate picture of the gross 
distribution. 

Temperature Sections (XBT) 

Horizontal sections of temperature at five selected depths from 
the XBT observations are plotted in figures 7a-e. Comparison of these 
to the previous plots demonstrates an obvious difference in detail between 
the results of the two surveys. In addition, the XBT temperature distributions 
appear to also give generally higher temperatures at a given location. 

Part of the reason for the greater detail in the XBT survey is 
obviously related to the increased number of observations, from 40 for the 


STD survey to 80. This increase greatly complicates interpolation during 
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plotting since one has too many points to obtain a simple mean distribution” 
but not enough to obtain a truly detailed distribution. Moreover, there is 
the limited accuracy of the XBT (? .20¢ and + 1% of depth) and the digitizing 
errors (+ .19C and * 1m) that accompany reading of the chart output. The 
net result is a number of closed contours which may or may not accurately 
represent the real oceanic situation. Although the latter case is most 
likely, the fact that the isolated contours exist mostly in the region where 
the eastward moving current divides (and is joined by the northward flow from 
the south along the coast -§ 2.6) does not completely rule out the existence 
of such closed isotherms at fixed levels. The comparatively large values of 
the XBT survey may also be related to its low accuracy. Other reasons that 
could account for the difference with the STD survey are: the difference 
of about | m between what is considered the 'surface' for XBT and what is 
considered the 'surface' for the STD; the tendency for the XBT thermistor 
to crack or be otherwise damaged and still produce what appears to be a 
reasonable trace. This 'recovery' of the XBT has the drawback that although 
the shape of the profile is real, an offset takes place whereby the absolute 
reading cannot be accurately determined. In this regard also, it should be 
mentioned that the known failure rate was about 20% with some XBT's showing 
malfunction immediately upon entering the water and others not until a 
few hundred meters. In the latter case, the temperatures to that depth 
were usually considered to be correct. 
It therefore seems reasonable to remark that an XBT survey is much 

more inaccurate than a corresponding STD survey in the sense of giving 

an ‘instantaneous' view of the spatial temperature distribution. There are 
just too many sources of error in this instrument to expect jt to give a 


detailed (< 0.20C) knowledge of the temperature structure, not withstanding 
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effects due to internal motions. 

Salinity Sections 

The low surface salinities appearing off the coast in figure 8a are 
indicative of the large amount of river runoff in this region that is occur- 
ring at this time of year. Included in this, of course, is that from the 
inland areas (via the Fraser River, for example) which is being advected 
seaward through Juan de Fuca Strait. At 50m, the salinity still decreased 
shoreward but with a much reduced gradient. By 100 m, the pattern had 
changed, with maximum salinities now appearing near the coast and Station P, 
and minimum values near 136°W longitude. At this depth, of course, the 
horizontal distribution is complicated the most by internal motions associa- 
ted with the halocline (~ 100-200 m) and undoubtedly accounts for the irre- 
gularity of the pattern. Below 100 m, salinity becomes increasingly uniform 
throughout the region although at 200 m more coastal values are Slightly 
larger than those in the interior (figure 8d). In addition, there appears 
to be a difference in the salinity distribution between the region west of 
128°W and that to the east, at least in the upper 100 m. Line P station 5 
(48°50'N, 128°40'W), therefore, seems to mark a transition from the coastal 
to the oceanic domain (see e.g. Dodimead, Favorite & Hirando, 1963). 

Density Sections 

Figures 9a-h show the horizontal distribution of sigma-t (s,) 
which is related to the density p at salinity = s, temperature = t and 
pressure = 0 db via the relation 


G, wu) lo 


£%.F Sue 9 ° 
As with the temperature and salinity from which it is derived, the vw, 
distribution has more or less a zonal orientation in the oceanic interior 


which gives way to coastally orientated contours in the near-shore region. 
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In all cases, there is a shoreward decrease in density, with maximum spatial 
gradients occurring essentially near the coast. 

This shoreward decrease in density clearly indicates that upwelling 
was not at yet taking place in the particular coastal region observed. 

At 100 m, the increased complexity of the Spatial pattern is most 
certainly the result of the irregular Salinity distribution Heese bed in the 
previous subsection. The strong dependence of density on salinity, therefore, 
also makes it impossible at this depth to obtain a consistent picture of the 
mean density field from only one set of observations. This dependence also 
accounts for Line P station 5 againmarking the transition from the coastal 
to the oceanic domain. 

Potential Energy Anomaly Sections 

In order to obtain a measure of the mean direction and magnitude 
of the baroclinic flow associated with the density structure, it is necessary 
to calculate certain properties which average over depth and include the 
influence of the earth's rotation. One of these is the 'dynamic height' aD 
defined at some pressure P by 

ADCP) = Js JP, 
where 6 is the anomaly of specific volume (density !). The spatial gradient 
of this quantity divided by the local value of the Coriolis parameter, f, 
then gives the magnitude of the surface geostrophic currents relative to 
the pressure surface P; the current direction is parallel to AD contours. 


In a similar manner, the potential energy anomaly y is defined by 
St x 
PRE) = g Jenene (1) 
r@) 


where g is the acceleration of gravity. Since 


"Neely Fag) of tA Denia 5 
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the spatial gradient of x divided by f gives the magnitude of total mass 
tramSport per unit horizontal width from the aufatt to the depth (z) at P. 
That is, relative to z, the vertically integrated geostrophic transport 
component U parallel to the x-contours is given by 

BAG 9e S315 te ne 
where s iS normal to the contours. 

Figures 10a-g are plots of ys width of 1 cm for some selected 
levels. Moreover, in the calculations of (1), the pressure has been replaced 
with 

P= pqz 
by means of the hydrostatic equation since the error in going from P to z 
is only about 1% at 1500 m. If we now assume that baroclinic motions at 
1200 m are very small, figure 10e can be used to determine the mean mass 
transport within the top 1200 m in the vicinity of Line P. Basically this 
transport approaches the coast from the west as a zonal flow with the more 
southerly portion eventually turning southward and the more northerly 
portion moving northward, in agreement with previous findings in this region 
(c.f. Dodimead, Favorite & Hirarmo, 1963). East of 130°W longitude, the 
transport apparently moves northward, parallel to the coast and with a much 
greater magnitude than in the more interior regions. The closed contours, 
if real, could be a consequence of the large amount of variability associated 
with the north-south splitting of the currents (Thomson, 1971). 

The magnitude of the transport between any two contours is readily 
Obtainable since it is simply the difference divided by the mean of the 
Coriolis parameter over the line of separation. At 50°N, for example, 
contours whose difference is 0.5 (x 108 ergs/cm) have between them a net 


is 


transport of 0.5 x 10 ~ gm/s_ or 0.5 sverdrups, referred to the appropriate 
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depth. The average speed of the surface current in cm/s’ relative to the 

depth of the potential energy surfaces can then be found by dividing the 
transport magnitude by the depth and spatial separation of the contours (in cm).— 
For contours separated by 25 km at 1200 m, mean surface speeds are about 

2 cm/s, while near Vancouver Island where the distance between contours is 

less, speeds increase to the order of 4 cm/s. 

Except for a slight decrease in magnitude, the mean transport 
relative to 1000 m is similar to that relative to 1200 m. At a reference 
level of 500 m there remains a good similarity with the deeper patterns 
although the magnitudes are much reduced. Finally, even referred to 100 m 
the transport retains the general pattern for deeper waters although it is 


relatively very small. 
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2.3. A surface temperature comparison 

The solid line in figure 11 is a plot of the surface (0 m) 
temperature obtained from the STD versus time of observation. Corresponding 
temperature values from XBT traces (x) and bucket samples (0) have also been 
included. 

Up to and including consecutive station 14, which marked commence- 
ment of the shore-bound portion of the cruise, the STD values are consistently 
lower than those obtained by the other two methods. During the inbound legs. 
however, the difference between the STD and bucket values appeared to fluc- 
tuate in sign as well as magnitude (no on-station XBT's were taken). These 
remarks are clearly demonstrated by figure 12a which gives the bucket 
temperature minus the STD temperature and the bucket temperature minus the 
XBT temperature versus time of observation. Furthermore, the curves show 
that the difference between the bucket temperature and the XBT temperature 
is consistently less than that between the STD and the bucket but that the 
fluctuations in the two cases remain in phase. 

Since May is a period of net solar flux into the upper layers of 
the subarctic Pacific, the above differences are for the most part related 
to the existence of temperature gradients in the upper 1 meter of water. 
Thus, the STD, which must be lowered at least a meter into the water to 
avoid spiking of the signal due to surfacing during ship roll, will measure 
a slightly different temperature than that from the bucket sample which 
essentially gives an average temperature in the top 30 cm. In the case of 
the XBT, its thermistor requires a finite time to respond to the ambient 
temperature and so can be expected to give a slightly different reading 
than either the STD or bucket sample. Nonetheless, figure 12a indicates that 


this response time is at least rapid enough to attain ambient temperatures 
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before a depth comparable to that of the 'surface' STD is reached. In fact, 
in five of the thirteen available cases the bucket and XBT temperatures were 
identical to within + 0.05°C which is rather surprising when ore considers that 
the individual accuracies of these instruments is at best t 0.1°C. 

On first considerations, the diurnal variation in the solar heating 
would be expected to account for a large portion of the difference in 
temperature measured by the three methods. Based on this effect alone, 
greatest differences should occur sometime in the afternoon when the surface 
layer is being heated the most. Contrarily, smallest differences should 
occur at night when radiative cooling is at a maximum. According to the 
LAYMORE data, this would appear to be the situation particularly from May 14 
Onward where positive peaks Occur regularly between noon and 1700 hrs_ PDST 
(= GMT - 7 hrs ). Similar peaks are also seen to occur in the air temperature 
at Bbut 10 m (figure 12b). Moreover, the negative values suggest significant 
cooling of the top thin layer usually takes place at night. The difficulty 
in relating cooling to a radiative effect, however, is that it is much too 
large in relation to the air-sea temperature difference -- the latter differen- 
ce being either small or negative (figure 13a). Also, at times prior to 
May 14, the succession of maxima and minima in figure 12a is far less 
regular than would be expected. The large air temperature of May 13, for 
example, iS associated with a negative bucket - STD value while the large 
positive peak on May 10 is not associated with a corresponding air temperature 
maxima. 

Some of the latter features may be explained if we now include the 
effect of the wind. The negative temperature differences in figure 12a 
which existed after May 15 are probably the result of the night-time surface 


cooling and mixing of the winds that had attained speeds greater than 10 knots 
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on these occasions (figure 13b). The relatively low daytime values during 
the 10th, 11th and 15th appear also to be related to wind strength maxima. 
The minimal temperature differences at 0700 GMT on May 9 and 2300 GMT on 
May 13, on the other hand, are not relatable to wind-strength but to the 
fact that the sky was completely overcast and a slight drizzle was falling -- 
the only two days on which the situation occurred. 
In connection with the wind-effect, it should be noted that in- 
creased winds also result in increased ship drift. Since all sampling 
was done on the weather side (starboard), there was, therefore, an increased 
amount of turbulent wake activity and hence mixing at the surface. This, 
together with the breaking of waves on the hull, would lead to measured 
surface values indistinguishable from those produced by the winds in the 
absence of ship interference. It should also be appreciated that instrument 
and digitizing errors could account for some of the observed differences. 
Finally, a comparison of humidity (figure 13c) obtained from the 
air temperatures of figure 12b with use of figure 12a, indicates a fairly 
good correspondence between the bucket - STD temperature difference and 
water vapor content. In particular, the trends in the two graphs are 
almost identical as might be expected for phenomena arising as a result 


of surface heating. 
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2.4, A Surface Salinity Comparison 
A plot of the surface salinity obtained from the STD and bucket 
samples is illustrated in figure 14a. The low values, including those 
centered around the 16th of May, are associated with the coastal zone off 
the B.C. coast and are a result of river discharge together with local 
advection. The effect of these is further indicated by figure 1] which shows 
that temperature maxima are coincident with each of the low salinity regions. 
For the most part, figure 14b indicates that the bucket salinities 
are larger than the STD salinities although differences less than 0.05°/00 
may not be significant because of the limited STD accuracy. Furthermore, in 
the near-surface layer of water, temperature gradients produce spikes which 
must be smoothed out. This introduces an additional uncertainty of about 
+ 0.01°/o0. Nonetheless, it would appear that the surface water sampled by 
the bucket was slightly more salty than that just below the surface. This, 
together with the absence of discernable transient haloclines in the top 
few meters, suggests that prolonged periods of evaporation which would account 
for such a situation were only beginning to take place. The large peak that 
occurs on the 17th of May, however, may be an example of persistent localized 
evaporation since it also corresponds to the occurrence of a positive | 


temperature gradient at the ocean surface (figure 12a). 
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2.5. Tidal/Inertial Motions in the Vicinity of Line P 

The subarctic Pacific water is typified by a large permanent halocline 
between about 100-200m over which the salinity increases from approximately 
32.5°/00 to approximately 34.0°/00 As a result, the 33.0°/o0 and 33.5°/00 
isohalines are located well within the depths of most vertically stable water. 
The time variations in the depths of these surfaces should, therefore, be 
indicative of any vertical motions occurring at their respective levels. 

Plotted in figure 15a are the depths of the above isohalines versus 
the observation time for the LAYMORE STD casts (those for the VANCOUVER and 
QUADRA are also included). Some of the variability is, of course, attributable 
to the spatial trend in the mean depth that was associated with all oceano- 
graphic properties in the region. Their relative shallowness in the coastal 
region, for example, was clearly contrasted to the deeper levels in the more 
interior region. In addition, however, there appeared to be a distinctive 
periodicity to the depth changes at periods one would expect for the tides 
(~12 and 24 hours) and inertial motions (~16 hours). Internal waves with 
periods of minutes to hours and wavelengths of kilometers would necessarily 
be superimposed on these longer period motions but could not be expected to 
produce such regular variations between stations which are separated by 2-4 
hours and 20-40 kilometers. On the contrary, such spacing would have been 
sufficient to obtain a nearly regular sampling of any tidal/inertial motions 
provided their amplitude exceeded those of the internal waves and the spatial 
trend in the depths of surfaces of constant property. 

The relative brevity of the record (~ 10 days) obviously limits 
the number of tidal constituents that can be identified. Moreover, the 
fact that the observations were not taken at a fixed location would cause a 


spatial aliasing of the temporal variations, which in turn would be dependent 
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on the direction of the ship's course in relation to the direction of tidal/ 
inertial wave propagation. Availability of computer programs to determine the 
spectrum and the harmonic tidal constituents, however, made analysis worth- 
while despite the limitations of the data -- the only step required to meet 
the input formats being a smoothing of the records to give values every 
hour (see figures 15b, c). | 

Figure 15d gives the spectra obtained for the unfiltered time vari- 
ations in depth of the 33.0°/00 and 33.5°/00 isohalines (a Fast Fourier 
Transform method was used to determine the necessary Fourier components). 
Aside from the large amount of 'energy' associated with spatial trends, there 


2 6 x 1072 and 8 x 1072 


are three main peaks near 4 x 10° cycles/hour corres- 
ponding to the diurnal, inertial and semi-diurnal periods. That the magnitude 
at the diurnal period appears to be greater than that at the semi-diurnal 
period may have been due to the separation of hours and tens of kilometers 
between stations whereby the spectra are determined by the sampling as well 

as by the actual tidal motions. 

The magnitudes of the various tidal constituents for the two iso- 
haline variations is presented in Table II (the inertial period is not included). 
Analysis has been based on the commonly used relation (Schureman, 1958): 

Gb = ee Pe cos tut tay 


© 


where h = depths of the isohaline surface in meters at time t, 
Ho = the mean depth over the record length, 
F = factor for reducing the mean amplitude H to the year of calculation, 
w = frequency of the constituent, 


t 


time reckoned from the initial epoch (midnight 1973), and 
Pe U = value of the equilibrium argument of the constituent when 


t = 0 (longitude dependent). 
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It is clear from this table that the Lunisolar tide K1 (1 cycle/day) 
and the Lunar tide M2 (~ 2 cycles/day) were the major constituents, with 
aliasing probably accounting for a large portion of the magnitude at the 
higher frequencies. The fact that the semi-diurnal component is of equal, or 
greater, magnitude than the diurnal component in this analysis (the reverse of 
the spectral analysis case) is most likely a consequence of the inertial 


variations having a period closer to 12 hours than to 24 hours. 
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2.6. Temperature Microstructure in the Vicinity of Line P: 
Evidence of the mixing of two water masses. 

Discussion 

At the time of the LAYMORE cruise, it became increasingly apparent 
_ that, for the deep (300-1500 m) section of the STD casts especially, there 
was a marked decrease in the amount of temperature fine-structure seaward of 
the near-coastal regions. This is examplified by figure 16 in which the deep 
cast at consecutive station 11 (49°41'N, 140°37'W) is plotted together with 
the deep cast at consecutive station 5 (48951'N, 128°940'W). The noisiness 
of the signal at the more coastal station was definitely much greater than 
that at the more oceanic station, with significant temperature structure 
still prevalent at 700 m in the former. 

The above discussion is of course only qualitative. In order to 
obtain a consistent measure of the amount of temperature fine-structure, it 
is necessary to do some sort of numerical analysis on the digitized records. 
One of the simplest and most straightforward forms of analysis 7s to determine 
the variance by first calculating a running mean and then subtracting it from 
the actual record. To this end, both the shallow (0-300 m) and deep (300- 
1500 m) parts of the casts have been smoothed using a box-car running mean 
of 19 m width in the upper layer and 95 m width in the lower layer. Choice 
of these averaging intervals was made after a number of trials to determine 
which gave the best visual fit to the actual temperature traces. As it 
turned out, an interval of 19 points gave the most satisfactory results, and 
since vertical resolution of the analog traces is 1 m and 5 m in the respec- 
tive layers, the above averaging depths were obtained. Subtraction of the 
mean values THs at each digitized point j, from the actual values Tj then 
gave the fluctuations T'; at each point (see e.g. figures 17a, 4). The 


variance 0% (where o is the standard deviation)over any block k was then 
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FIGURE 16: "LAYMORE"S.T.D. profiles for two Line P stations 
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determined via the relation, 
a = ct at ~ 

where N is the number of data points within the block. In these calculations, 

| a block consists of values from a specified depth range for each station. 

Upper and lower layer blocks 

For these cases, the variance has been determined for the two depth 
intervals 0-300 m and 300-1500 m corresponding to the two ranges used during 
lowering of the STD. The calculated values in °C2 are plotted in figures 
18a, b with the numbers in brackets representing results from the second 
(later) occupation of a given Station. In the top 300 m, there was obviously 
a sharp increase of the amount of temperature structure in the outbound 
compared to the inbound observations in the coastal region. Within this 
region, therefore, inbound values rather than outbound values have been used 
Since in time they correspond more closely with those north and south of the 
Line. West of station 6 either value isapplicable since the contour positions 
are nearly the same for both. As figure 18a shows, the contours indicate a 
maximum region extending from about 47°N, 130°w through to 49°N, 128°W with 
an additional relative maximum occurring 2° south of station 4. Relatively 
large values also occur west of the northern end of Vancouver Island to about 
140°W longitude. 

The difficulty with associating these results with temperature 
‘microstructure', however, is that the smoothing process is incapable of 
removing all temperature variations with depth which are due to large scale 
features (10's of meters), such as the surface mixed layer. Therefore what 
one has is a combined measure of the detailed fine-structure and the magnitude 


of the large scale temperature gradients. Moreover, the effect of the latter 
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features can be highly variable over short periods because of internal waves 
whose motion can increase or decrease the local temperature gradients by 
compressing or expanding the depth between isotherms. In the latter situation, 
application of the low pass running mean will give rise to larger apparent 
fluctuations than in the former situation. Nonetheless, the present plot does 
give an indication of the overall temperature structure and confirms the 
notion of greater temperature variability in the coastal areas relative to 
the deeper oceanic areas. Furthermore, all indications of the amount of 
microstructure in the top layer are not completely lost since by dividing this 
layer into two, the large scale temperature features become confined to the 
upper 150 m only. This will be done shortly. 

With regards the temperature variance in the 300-1500 m range, one 
can be more assured that values obtained are representative of the fine 
structure alone since the mean profile remains smooth throughout. This is 
supported by figure 18b which indicates relatively small changes in variance 
between temporally spaced observations at the same location. For this range, 
therefore, there is little difference in contour positions if the inbound 
casts rather than theoutbound casts are used for determination of the variance. 
Comparison of figures 18a and b shows that the spatial distribution in the 
deeper range is nearly identical to that in the 0-300 m range although the 
‘magnitude is reduced by an order of ten or so at all positions. Thus, in 
this layer we again find maximum variances occurring in a narrow strip near 
130°W, with values west of this decreasing monotonically towards Ocean 
Station P. It should also be noted that the value of 1.01 used for station 5 
probably gives an underestimate of the values in this region since it is from 
the outbound portion of the cruise while all others are from the inbound 


portion. Judging by the increases that occurred at stations 4 and 6 a value 
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closer to 1.5 might be more realistic. A similar situation arose in connection 
with the 0-300 m range, but there it was expedient to simply ignore the out- 
bound value at station 5. 

Aside from what occurred in the top 150 m, which will be discussed 
| in the next subsection, the most plausible explanation for the proceeding 
distributions is that the positions of large variances are associated with areas 
of horizontal mixing between two water masses of different temperature. The 
larger variances mark the 'front' where the confluence of the two waters is 
taking place. According to the potential energy anomaly plots of figures 
10a-q and the T-S diagram of figure 19, the mixing is in fact between the 
relatively cool mass of water moving eastward from the oceanic interior and 
the relatively warm mass of water moving northward off the coast of BC. and 
Washington. 

Division into nine blocks 

In order to determine the degree of temperature variance with 
depth, the two previous ranges have been further divided into the following 
nine layers: 1) 0-50 m, 2) 50-150 m, 3) 150-300 m, 4) 300-500 m, 5) 500-700 m, 
6) 700-900 m, 7) 900-1100 m, 8) 1100-1300 m, and 9) 1300-1500 m. Layers 1 
and 2 have been so chosen in an attempt to separately confine the effects 
of averaging through the mixed layer (< 50 m) and the dicothermal layer 
( 50-150 m), respectively. Figure 20a gives a measure of the magnitude of 
the temperature gradients associated with the surface mixed layer. Obviously, 
the gradients are much larger in the coastal areas than in the oceanic areas. 
There is also a greater amount of temporal variability in the temperature 
Structure in the coastal region. For the 50-150 m range, the temperature 
variance is not so dependent on relative proximity to the coast although, 


for the available stations, the greatest time variability appears to have 
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occurred in the coastal region. Consideration of the original data suggests 
that the latter result is probably indicative of steepening of the vertical 
temperature gradients by surface mixing and internal waves at places where 
there was much temperature structure to begin with. 

The temperature variations plotted in figure. 20c for the range 
150-300 m conform more to those for the deeper ranges than for the more 
shallow ones. Therefore, for depths exceeding 150 m the values obtained 
for the variances are mostly independent of contamination from large scale 
gradients that are not removed by the smoothing technique. What remains 
is microstructure produced by the confluence of the two water masses. 
Figures 20d-i then give the following picture for 150-1500 m range. 

a) in the range 150-700 m, maximum variation occurs near the coast 
in a region extending from about 47°N, 131°W to about 49°N, 129°W. Westward 
of the coast, variances tend to decrease, with minimum values appearing in 
the vicinity of station 12 (Line P). In addition, the horizontal sections 
indicate a fairly extensive area of relatively large values near 500N between 
Vancouver Island and about 135°0W. Vertically, the maximum values for stations 
east of the longitude of Line P station 7 occur always in range 4, 300-500 m 
(see figure 2la). West of this (including station 7), however, maximum 
values with depth alternate between levels 3,4 and 5, with the exception of 
those at the longitude of Line P station 8 which again have their maximums 
at level 4 only (figures 21b, c). 

b) in the range 700-1300 m, the region of maximum structure has 
shifted more towards the coast although values to the north of 50° continue 
to remain relatively large for these depths. In the vertical, the magnitude 
of the variance is nearly uniform with depth at all stations and is small 


(<.5.x 10720¢2), 
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c) averaged from 1300-1500 m, the magnitude of the temperature 
variance is comparable to that in (b) but with the maximum in the horizontal 
shifted more away from the coast. The pattern is similar to that in (a). 

Because of the continued persistence of the strip of maximum 
variance near 130°W longitude in all horizontal sections, it would appear 
that the northward flowing coastal current extends below 1500 m. The fact 
that the magnitude of the variance generally decreases with depth is simply 
due to the corresponding decrease in the difference in mean temperature 
between the two mixing water masses (see figure 19). The most obvious 
exception to the latter trend is the finding that below 150 m the greatest 
variances east of 132°W longitude occur between 300-500 m. This Suggests that 
the temperature contrast between the two water masses was greatest at this 
level and/or that the mixing was particularly vigorous. Further seaward, 
where the amount of temperature structure with depth is significantly 
smaller, the level of maximum variance is much more variable. Such a result 
is to be expected for this oceanic region since its temperature structure 
is an inherent feature of a particular water type formed on the western side 
of the subarctic Pacific and not a result of more local mixing between two 
water masses. 

Another feature of the variance distribution is that, compared to 
other levels, there is in the general range 700-1300 m a greater westward 
penetration of the northward flow with maximum extent occurring near 48°N 
latitude. Below 1300 m, however, the pattern outlining the mixing region 
returns to the form for depths less than 700 m. Associated with this 
changing spatial distribution with depth is the varying magnitude of the 
temperature variance at consecutive station 31 (47°10'N, 130940'W). Exclud- 


ing level 2, the variances at this station for depths less than 700 m appear 
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to be much larger than at any neighbouring station and usually represent the 
maximum value for the survey. It is, therefore, not unlikely that these 
rather large values are somehow related to the presence of Cobb Seamount 

(~ 46945'N, 130940'W) only 35 km to the south. In particular, there may be 
some sort of turbulent wake associated with this mountainous feature since, 
by rising to within 35 m of the surface from depths of over 2700 mina 
distance of only 20 km (see figure 22), it penetrates abruptly into most 
levels of the baroclinic flow. In the presence of such a wake, the horizontal 
mixing of the two different temperature water masses in this region may be 
enhanced compared to more distant regions. This would be particularly true 
of the nove Shallow flow where the baroclinic motions are strongest. A 

less speculative understanding of the influence of the seamount, however, 


awaits a more detailed survey. 
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FIGURE 2la: Temperature variance with depth: Consecutive stations 
4(e), 5(¢), 34(4), 36(x ): 
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FIGURE 2lb: Temperature variance with depth consecutive stations 
6 (e), 27 (4); 31(x); 32 (¢) 
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3. The Variability of Water Properties Along Line P 

3.1. General discussion 

In this section, results of observation from the LAYMORE are used 

in conjunction with those from the weatherships VANCOUVER and QUADRA to 
determine the variability of water properties along Line P from 8-18 May. 
A plot of station number and longitude versus time for each of the three 
ships was presented previously in figure 2. One interesting feature of this 
was that the mean outbound speed of the LAYMORE was only slightly less than 
that of the VANCOUVER. Therefore, the relatively slow between-station speed 
of the former, which is only 11 knots compared to about 16 knots for the 
weatherships, was nearly compensated for by the shorter time required to 
do each station. As indicated in the Introduction, this is not a reflection 
on the ability of the crew or scientists but rather on the ship designs 


which make observations from the LAYMORE easier to perform. 
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3.2. A Comparison of Surface Salinities and Temperatures 

The surface salinities and temperatures from the STD observations 
of the three ships is shown in figures 23a and 23b, respectively. Station 5 
is seen to delineate the transition region from the more uniform structure 
in the oceanic domain to the more rapidly varying structure in the coastal 
domain. This is mostly true of salinity, where there is a rapid salinity 
increase seaward of the coast and a nearly uniform and constant salinity 
in the interior. Because of the small changes in the latter region, it is 
not possible to distinguish actual variations in space and time from those 
due to instrument error. East of station 5, however. «the variations are 
probably significant. With the VANCOUVER observations (x) taken consistently 
between those of the LAYMORE (.) and QUADRA (0), these variations indicated 
a shifting of the isohalines by horizontal currents and/or by changes in 
the amount of river runoff. 

Aside from instrument errors, the most likely explanation for the 
generally increased surface temperature over the 2-week period is that the 
LAYMORE observations were taken during passage of a fairly intense storm 
front. Asa result of cooling and wind mixing at that time, surface tempe- 
ratures taken by the LAYMORE would be somewhat low compared to those of the 
QUADRA and VANCOUVER which were taken-at more calm periods. This would also 
account for the greater temperature changes that occurred west of station 7 
as compared to more easterly stations since it was there that the maximum 
effect of the storm was first encountered (see also figures 25a-q). Northward 
advection of relatively warm water could account for some of the general 
temperature increase but not all of it since surface flow tends to be zonal 


and weak. 
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3.3. Depth of the Surface Mixed Layer 

During the spring, the accumulated effects of surface heating and 
wind mixing over most of the subarctic Pacific combine to produce a surface 
layer of water (< 50 m) that has a nearly uniform temperature with depth 
and which is truncated by a temperature discontinuity. Formation of this 
layer is part of a cycle that consists of an upper layer of uniform temperature 
to about 100 m in mid-winter (January), to highly temperature-structured upper 
waters in Summer when surface heating is greatest and wind least. Since in 
May the mixed layer depth is intermediate between these two extremes, it 
should be highly succeptible to variations induced by the winds through 
turbulent mixing. Moreover, as the pycnocline (which is essentially salinity 
determined) is always deeper than the mixed layer, the effectiveness of the 
mixing is not dampened by large scale density gradients. Hence, the mecha- 
nical energy input by the wind should be rapidly distributed to the base of 
the layer by turbulent diffusion where it will work against the temperature 
gradient. As a result, it should be possible to associate large temporal 
changes in mixed layer depth to the wind strength. Internal waves, tides, 
inertial waves and Ekman divergence (convergence) generated by the surface 
wind-stress curl will, of course, also influence the observations although 
only the latter is actually capable of altering the 'mean' depth in any 
spatially consistent manner over periods of days. According to Denman (1972), 
however, who conducted studies at Station P, the latter is probably not 
significant when compared to the turbulent mixing. 

Illustrated in figure 24 is the spatial variation in mixed-layer 
depth as measured upward from the top of the thermocline for the three 
sampiings of Line P. Between stations 1 and 7, the LAYMORE depths were 


consistently much more shallow. than those from the QUADRA, while those 
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from the VANCOUVER were intermediate. West of station 7, LAYMORE values 
deepened, while the weathership values rapidly sha] owed) A thin mixed 
layer was then recorded by the LAYMORE at station 12 and a deep mixed layer 
by the VANCOUVER at station P. As the dates on the curves indicate, LAYMORE 
' observations proceeded all others in every case while VANCOUVER observations 
proceeded those of the QUADRA except at stations 12 and P (see also figure 2). 

The surface pressure, winds, cloud cover and frontal positions at 
midnight (OOhrs) GMT for the period May 2 to May 18 are presented in 
figures 25a-q. On May 2, the cyclone which had started to form in the vicinity 
of the Aleutian Islands at the end of April, began to move eastward over the 
northeastern Pacific. East of about 150° longitude, however, winds remained 
light and skies clear. By May 5, the storm had penetrated only to about 
137°W .and was much reduced in strength while on the next day, its effect was 
limited to the vicinity of the cold front marked in figure 25e. By May 7, 
winds over the region of Line P were again light and skies partly clear. 
This situation prevailed until May 9 when “dicaeWen eastward moving disturbance 
had begun to penetrate east of station P. 

The fairly extensive period of calm weather that existed over 
Line P up to May 9 would then account for the shallowness of the mixed layer 
observed by the LAYMORE for stations 1 to 4. On the other hand, the slight 
deepening that occurred at eiesiors 5-7 appears to coincide with the time of 
arrival of the storm's edge at these positions. By the time the LAYMORE 
had arrived at station 8, the mechanical wind-mixing effect of the storm 
had had sufficient time to markedly deepen the layer. This would account 
for the abrupt deepening observed between stations 7 and 8. The shallowness 
of the mixed layer at station 12 (OOGMT May 12) resulted from the low winds 
(< 10 knots) and clear skies in its vicinity for the day or so prior to 


observation. 
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By May 13, the greatest effect of this particular storm was being 
felt east of 135°W on Line P and would account for the rapid deepening of the 
mixed layer observed by the outbound VANCOUVER. West of 135°0W, there was not 
a Significant increase in the wind until May 18 when the third storm was just 
beginning to move past Station P. This would explain the relatively shallow 
values from station 8 westward during the VANCOUVER cruise as well as the 
large depth observed at P at 1820 on May 18. Finally, the continued deepening 
of the mixed layer in the more coastal regions as obtained by the QUADRA, 
are probably a result of the low pressure regions that built up near the 
coast after May 16. 

In general, then, it appears that the depth of the mixed layer over 
the period of observation is closely related to the availability of mechanical 
energy from the winds. Moreover, Ekman layer divergence (convergence) gene- 
rated by a positive (negative) wind-stress curl at the surface is clearly 
not of importance since it would produce an effect opposite to that observed. 
For example, cyclonic winds which would generate a divergence of the surface 
waters and therefore an upward vertical velocity at the base of the main 
Pycnocline should be associated with a decrease in the thickness of the 
upper layer; anticyclonic winds on the contrary should be associated with 
an increase in layer thickness. That the reverse Situation actually occurs 
is evidence of the unimportance of the Ekman divergent effect compared to 
the mechanical mixing of the surface winds in altering the thickness of 
the surface mixed layer. Below the depth of direct surface influence, 
however, the vertical velocities induced by the Ekman layer divergences can 


be of importance as shown in the next subsection. 
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3.4. Depths of the 33.0% and 33.5% Salinity Surfaces Versus Longitude 
The depths of these two surfaces as determined from the three 
samplings of Line P are plotted versus distance in figures 26a, b. Although 
much of the station-to-station change within a given cruise is attributable 
to internal motions (see § 2.5), it is apparent that there was an overall 
increase in the mean levels between the LAYMORE observations (9-12 May) and 


the weathership observations that followed. There did not, however, appear 


to be any consistent change in the mean level between weathership cruises. (Fig. 


The spatial trend of shallow values near the coast tending to 
deeper values seaward is typical of Line P and is simply a manifestation 
of the shallowing topography and coastal currents. The temporal increases 
in depth, however, must be the resultant of large scale processes over this 
corner of the Pacific. One possible explanation is that between the time of 
the LAYMORE and the weathership samplings there was a general relaxation in 
the upward vertical velocities at the base of the halocline as a result of 
decreased cyclogenesis over the region. Examinatton of figures 25a-q indicate 
that this may indeed have been the case since the cyclonic wind patterns that 
existed over much of the subarctic region from 9-12 May gave way to anti- 


cyclonic winds from 12-15 May. The resultant “spin-down" associated with 


26c), 


the Ekman layer convergence, which occurs within a half-pendulum day (~ 18 hours), 


would then allow sinking of the surfaces of constant physical properties to a 
new equilibrium level. 

Another feature of the subarctic Pacific which would induce similar 
results, would be a weakening of the Alaska gyre by the advent of anticyclo- 
genesis over the Gulf of Alaska. Its effect would be contrary to what is ob- 


served, however, since it would be expected to produce a relative decrease in 


the depth of the salinity surfaces along the Line through a decrease in northward 


advection into the gyre in a region where the depth of salinity surfaces increases 


southward. 
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3.5. Depth of the Dichothermal Layer 

As defined by Uda (1963) theDichothermal layer is a tongue of 
relatively cold water which sometimes exists at a depth of around 100 m in 
the subarctic Pacific Ocean (figure 27). On the eastern side of the ocean, 
Feeian of this water begins with intense winter cooling and wind mixing. 
By late winter this results in a surface layer (~ 100 m) of isothermal 
water that is colder than that below. With the advent of spring, surface 
heating occurs and warmer water is mixed downward but not far enough to remove 
the positive temperature gradient at depth. The feature then usually persists 
until the next cooling period (fall-winter) when it is either destroyed or 
enhanced. 

The difficulty with determining a depth for the dichothermal layer 
is that various internal wave modes can greatly distort the temperature 
profile. This distortion is particularly prohibitive in 0-300 m STD casts 
because of the expanded scale so where possible 0-1500 m reverse traces have 
been used. Even then, errors of order * 10 m (or more) in reading. the mean 
depth are readily possible. Furthermore, what can be interpreted as a 
didiothermal layer may in fact be associated with the mixing of different 
water masses and not with the past effects of winter cooling. As discussed 
in 2.6,this is particularly true for near coastal stations and depths 
greater than 150 m. Despite these, and the internal effects of vertical 
displacements, figure 28 gives a reasonably compatible set of plots of the 
depth of adichothermal-like layer versus station position. The shoreward 
shallowing of the layer is associated with a similar decrease in the depth 
of the top of the pycnocline since, as expected for a surface influenced 
feature, the depth of penetration of the winter overtuning and cooling is 


clearly limited by the stability of the water column. Again, station 5 
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appears to mark the transition from the coastal to the oceanic domain. In 

the latter region, the maximum depth of late winter cooling is seen to be 

around 100-125 m in agreement with Uda. 

The temperature gradient of thedichothermal layer defined by 
Bh Oey indie’ OZ os 

where ST is the temperature difference and Az is the thickness,-is given in 

figure 29 for the three sets of observations. Irrespective of the large 

amount of variation at each position due to stretching and compressing of 

the isotherms by internal wave motion, there appears to be a definite maxi- 

mizing of the gradient near stations 3 and 4. ° Categorically, it is not 

possible to state whether this jis a result of water mass mixing as described 

in §2.6. or whether it is the result of particularly vigorous mixing that 

occurred in the winter of 1971-72. Possibly it is a combination of both, 

although the fact that this region compared to others has relatively large 

temperature gradients to 1500 m, suggests the former. Therefore, the plots 

can be considered to identify the depth and temperature gradients of the 
dichothermal layer except in the vicinity of the coastal stations where 

the feature is masked by the horizontal confluence of the coastal and interior 


water masses. 
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3.6. Geostrophic Mass Transport Across Line P 

The concept of the geostrophic transport between two stations being 
proportional to their difference in the potential energy anomaly is based on 
steady-state calculations. Nonetheless, it is possible to obtain some idea 
| of the flow variability with time if the adjustment period of the density 
structure is short compared to that between observations. 

As indicated previously by the potential energy plots of figures 10a-9, 
the water motion in the vicinity of Line P would not be expected to have a 
constant component perpendicular to the Line. This is supported by figure 30 
which gives the spatial variation of the potential energy relative to 1200 m 
for each of the three transects (in each case, the component of integrated 
transport at right angles to the Line is such that the large values are to the 
right facing in the direction of the flow). Subtracting the value at 
station 12 from that at station 4 for each cruise further indicates a variation 
in the mean total northward transport from 5.76 sv (sverdrups) during the 
period of the LAYMORE crossing, which decreases to 4.77 sv during that of 
the VANCOUVER, to 6.57 sv during the QUADRA cruise. The greatest amount | 
of temporal variation in any one locality is between 128° and 135°W longitude 
where the eastward moving waters are being most influenced by the coastal 
regime. In the latter, between-station transport is maximum and directed 
parallel to the coast. Finally, taking 5.5 sv as the mean mass transport 
between stations 4 and 12 gives an average surface current speed v of 

V~ 0.43 cm/s 

or about 0.4 km/day at right angles to the Line. This compares to a mean 
northward speed of | 


Vv ~ 62.63 cm/S 


between stations 4 and 5 alone, or a surface drift of about 2.3 km/day. 
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3.7. Comparison of Temperature Variance Along Line P 

The temperature variances versus longitude along Line P for the 
three separate cruises are plotted in figures 3la-g. The first two give the 
block averages for the ranges 0-300 m and 300-1500 m while the remaining 
plots give the variances in five two-hundred meter ranges from 300-1300 m. 
The last range 1300-1500 m is omitted since it is very similar to the 
1100-1300 m Ss 

Averaged over the top 300 m, there is again the tendency for the 
near coastal values to be larger than the more oceanic regions. There is 
also a general increase in the coastal variances with time since chronolo- 
gically the order of ship sampling east of station 6 was LAYMORE (out), 
VANCOUVER, QUADRA and LAYMORE (in). This increase, moreover, corresponds 
almost exactly with the deepening of the mixed layer with time in the coastal 
areas (figure 24). The indication then, is that the apparent increase in the 
amount of temperature microstructure in the .0-300 m block is actually a 
manifestation of the increased large scale vertical temperature gradients 
rather than the result of a change in the horizontal mixing of water masses. 
Such is not the case for the 300-1500 m block. The significantly large tempe- 
rature variances east of Line P station 5(128940'W) together with the consis- 
tently smaller values to the west are a good indication of the persistence 
of the mixing process off B.C. The changes between cruises in the coastal 
region are then the result of spatial shifting of the regions of confluence 
and perhaps of the degree of mixing. 

Within each of the five two-hundred metre intervals, the spatial 
distribution of the variances closely emulates that of the average over 
these layers with maximum values during a particular cruise occurring in 


level 4 or level 5: The one exception is for the 500-700 m level(5) for 
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the QUADRA cruise which indicates relatively low coastal values and a marked 
maximal value at station 8 (134940'W). The latter is perhaps a manifestation 


of the discrete spatial sampling in which the region of maximum horizontal 


mixing was missed. 
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4, SUMMARY 

This study described some results from the data collected by three 
ships in the vicinity of Line P during May 1972. Although the oceanographic 
region and time covered were rather limited, a number of findings have been 
suena which have a much more extensive applicability in the subarctic 
Pacitic. 

Section 2 considered the data from the CFAV LAYMORE cruise alone, 
with the following results: 

1) Horizontal sections of temperature and salinity and hence 


density (*,) showed a marked zonal distribution of the isopleths in the 


+) 
interior oceanic region which became more orientated with the shore as the 
coastal domain was approached. 

2) The tendency for the contours of equal property to align them- 
selves with the coast was apparent at all depths in which there was a 
Spatially discernible difference in the magnitude of a given property. 

3) The transition from the coastal domain to the oceanic domain 
along Line P appeared to occur in the vicinity of Line P station 5 
(48051'N, 128°40'wW). 

4) In addition to the re-orientation of the contours as the 
coast was approached, there was a sharp increase in the spatial gradient and 
the time variability at most depths. 

5) According to the potential energy anomaly sections, the mean 
northward mass transport in the coastal domain was much more intense than 
in the more oceanic regions. Moreover, the commonly observed splitting of 
the eastward flowing West Wind Drift in the latter region into northward 
and southward tending flow seemed to be taking place to south of Line P. 


In particular, evidence of the splitting in the limited area covered by the 
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LAYMORE survey was confined to the region south of latitude 489N and west 
of longitude 134°W. 

6) Comparison of the STD temperature survey and the XBT survey 
indicated a significant difference in detail as well as in actual value at a 
given location. In part, this was a result of the twice as many XBT as STD 
observations. Mostly, however, the difference was due to the much greater 
inaccuracy of the XBT's compared to the STD and the tendency of the XBT's 
to malfunction but still give believable traces. 

7) Comparison of the surface temperatures obtained from the STDs 
XBT and bucket samples showed that, in most cases, the XBT values were closer 
to the bucket temperatures than were the STD values. This was associated 
with the fact that, in the surface layer in May, there is usually a decrease 
in temperature with depth. The STD, which must be lowered to about one meter 
to avoid surfacing of its sensors, then usually recorded a lower temperature 
than the other two. 

The exceptions, which occurred frequently during the LAYMORE cruise, 
appeared to be associated with times of reduced insolation (cloud cover; 
rain) and winds of 10 knots or greater, thereby implying an increase of 
temperature with depth in the surface layer at these times. 

8) The very small difference (£0.05°/00) between the STD salinity 
and the bucket sample salinity indicated then near absence of any salinity 
gradients with depth at the ocean surface. Moreover, this implies that, in 
general, very little evaporation or precipitation was taking place at this 
time, although the fact the bucket values were consistently larger than 
those of the STD may have been indicative of a slightly evaporative situa- 


ctor. 
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9) Time variations in the depths of the 33.0°/00 and 33.5°/00 
dsohalines surfaces have indicated that the accumulated vertical displacements 
of surfaces of constant property over periods of hours were tidal/inertial in 
nature. The Lunisolar tide Kl ( ~1 cycle/day) and the Lunar tide M2 
(~~ 2 eycreaydend were the major constituents, with some magnitude, at 
other frequencies due to aliasing as a result of the brevity of the record 
(~10 days). Maximum vertical displacements associated with the combined 
tides and inertial motion were approximately 25 m. 

10) In the upper 150 meters of water in the subarctic Pacific, 
the horizontal distribution of the variance associated with the temperature 
structure in the STD traces was closely related to the spatial variation in 
the intensity of the large scale temperature gradients (a manifestation of the 
low pass filtering technique). During the LAYMORE survey, for example, 
the spatial distribution of the temperature variance for the depth interval 
0-50 meters appeared to be a measuré of the spatial variation in magnitude 
of the thermocline at the bottom of the surface mixed layer. For the 
50-150 meters depth interval, it appeared to be a measure of the spatial 
variation of the temperature gradients associated with dicothermal layer 
(relative temperature minimum layer). 

In both ranges, greatest magnitude and temporal change occurred 
near or within the region of transition from the oceanic to the coastal 
domain. Minimum values and minimum time variation occurred in the oceanic 
domain. In the case of the 0-50 m layer these findings were indicative of 
the greater amount of wind-mixing that occurred in the coastal region compared 


to the oceanic region over the period of the LAYMORE cruise. 
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11) Below 150 meters, the temperature variance calculated from the 
vertical STD traces gave a more accurate measure of the amount of 'micro- 
Structure’ than for shallower depths. That is, the magnitude of the tempe- 
rature fluctuations about the mean at the former depths was determined by 
small-scale internal processes rather than by the large scale surface derived 
Processes and the spurious values resulting from the filtering. 

In particular, the horizontal distribution of the variance in the 
region studied, delineated the area of confluence between the warm northward 
moving water off the BC-Washington coast and the colder eastward moving water 
of the West Wind Drift. As shown, this mixing region was situated about 
200-300 km offshore. Between 150-700 m, it extended more or less northeast 
across the area near 48°-49°N and 130°W while between 700-1300 m it has shifted 
more towards the coast except for a tongue-like penetration to the west near 
480N. Despite the much reduced magnitude, the variance at 1300-1500 m 
continued to define the region of confluence but with a pattern similar to 
the former range (150-700 m) rather than latter range. A second geographic 
area of relative maximum temperature structure also appeared south of the 
Queen Charlottes at most depth levels. Seaward of these regions of confluence, 
the variance decreased rapidly to Significantly smaller values. Minimum 
variances occurred invariably in the vicinity of stations 12, 13 and 14 
(~ 143°W longitude). 

The much greater temporal variability in the variances in the 
coastal region than in the oceanic region was indicative of the shifting of 
the 'strip' of confluence and/or the intensity of intrusion of one water type 
into another off the coast. In the oceanic domain, the temporal changes in 
the variance are mostly associated with much smaller scale processes such as 


the breaking of short internal waves and cannot, therefore, begin to produce 
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temperature fluctuations of the same order as those due to the mixing of 
the two different water masses. 

12) Considered in the vertical, the variances associated with the 
temperature structure below 150 m indicated a maximum at level 4 (300-500 m) 
for the geographic area east of Line P station 6 (49°02'N, 130°40'W). West 
of station 6, the level of maximum variance varied between level 3 (150-300 m), 
level 4 (300-500 m) and level 5 (500-700 m). In the latter area, of course, 
maximum variances were much less than those in the more coastal areas except 
for the deepest levels where values became small everywhere. 

The consistent maximum at level 4 for those areas east of station 6 
is indicative of a greater temperature contrast between the water masses at 
this level or of the greater degree of mixing, or both. The spatial varia- 
bility in the depth of maximum variance for the areas west of station 6 
is indicative of the much smaller and more locally determined temperature 
structure found there. 

Section 3 considered the data from the two weathership surveys 
along Line P in May of 1972 together with that from the LAYMORE (outbound), 
with the following results: 

1) A comparison of surface STD salinities for the three separate 
samplings of Line P showed that differences above instrument error occurred 
in the coastal domain only. The changes that did occur in the latter 
region were most likely associated with the effects of river runoff and/or 
variations in the local current patterns. 

2) A comparison of surface STD temperatures indicated a general 
warming of the water with time along Line P. This ae particularly true 
of the period following the LAYMORE observations during which winds were 


generally lighter and insolation greater than previously. 
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3) The temporal variation of the depth of the mixed layer at a given 
location during May was directly associated with corresponding variations 
in the magnitude of the surface wind-stress. Moreover, the response time 
of this depth to the wind appeared to be less than one day. Vertical 
velocities associated with a divergent (or convergent) Ekman layer, if 
present, did not seem to be of any importance in altering the depth of the 
mixed layer. Therefore variations in the intensity of turbulent mechanical 
mixing of the winds and not vorticity effects were responsible for changing 
the layer thickness. 

In addition to the temporal variations, there was a distinctive 
shallowing of the mixed layer depth near the coast. This is related to the 
general upward displacement of all property surfaces in this region which, 
through increased water stability at shallower depths than in the interior 
regions, limits the depth of penetration of the mechanical wind mixing. 

4) Qualitatively, the overall increase in depth of the 33.0°/00 
and 33.5°/00 isohalines between the LAYMORE and weathership cruises along 
Line P appeared to have been related to a general anticyclogenesis in the 
wind activity over the region as a whole, rather than to the local wind 
magnitude. This would suggest that these isohalines, which lie well within 
the permanent halocline, have had their equilibrium level deepened through 
a decrease in the upward vertical velocity induced by the divergent surface 
Ekman layer usually associated with this part of the subarctic Pacific. 
Station-to-station fluctuations produced by internal waves (tides) were, of 
course, present but could not under most circumstances be responsible for the 
general shift of the mean levels. 

Between the weathership samplings of Line P, the situation was 


less well defined although, in general, the depth for the two isohaline 
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surfaces appeared to have been slightly greater during the VANCOUVER cruise 
than during the QUADRA cruise. 

5) The existence of adichothermal layer along Line P indicated 
that fairly intense surface cooling had taken place in the winter of 1971-72. 
The shoreward decrease in the depth of this layer was again associated with 
the upward displacement of all property surfaces near the coast since the 
depth of penetration of winter cooling and wind mixing would have been 
limited by the stable water of the halocline. 

Maximum temperature 'gradients' associated with the dichothermal 
layer occurred at Line P stations 3 and 4, which may be due to a greater 
intensity of winter cooling there or to the fact that these stations were 
Situated within the region of confluence between two water masses of different 
temperature. 

6) Comparison of the potential energy anomalies referred to 
1200 m along Line P for the three surveys, showed a variation of 5 to 63 
sverdrups in the total across-the-line mass transport. Maximum variability 
with time between the cruises occurred between 128°-135°W longitude in the 
region where there existed the greatest tendency of the West Wind Drift to 
be forced to flow parallel to the coast. 

7) Consistent with the results of the LAYMORE survey alone, 
weathership observations along Line P verified that the maximum variability 
in the STD temperature profiles existed in the vicinity of stations 4 and 5. 
Seaward of these stations there was a rapid decrease in the variance at all 
levels. As before, this implies the existence of mixing between a northward 
flowing coastal current and the eastward moving West Wind Drift off the coast. 
The significant variation with time of the coastal variances is a clear 


indication of the shifting position of confluence. 


Too 


8) The combined results of all three ship surveys further confirms 
the location of Line P station 5 as marking the transition from the coastal 


to the oceanic domain along the Line. 
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INTRODUCTION 


This summary of hydrographic and oceanographic information on. some British 
Columbia estuaries was initiated for the Pacific Region Board Working Group 
on Estuaries. It is concerned with those places on the Pacific Coast of 
Canada where industrial development and population growth are closely 
associated with river estuaries. The data sources include publications of 
the Fisheries Research Board of Canada, the Marine Sciences Directorate, 
Fisheries Service, and the University of British Columbia Institute of 
Oceanography. A great deal of the information was discovered upon personal 
enquiry to scientists and technicians in various government departments. 
Internal publications of these departments are often not catalogued in a 
fashion that makes them available to those people in other areas who could 
make use of them. A great quantity of data also remains on file, unpublished 
and thus not easily accessible. With changes in personnel and passage of 
time, data on file becomes less available and often later workers have an 
insufficient basis for evaluation of the information. 


Owing to the difficulty of searching many diverse sources and the impossi- 
bility of communicating with every person working on estuaries in the many 
concerned agencies, the present collection of information cannot be considered 
complete. Further contributions and comments from those more intimately 
involved in some particular estuary are solicited and will serve to bring 

the study closer to being a complete survey. 


The catalogue covers the estuaries of the rivers listed on page 5, in the 
order shown. Their location is shown on the first chart. For each estuary 
there is a chart with the area of concern outlined and the following informa- 
tion: 


Hydrography 


- Date of the most recent survey. 

- Scale of field sheets (F.S.) available. 

- Chart numbers and scales. 

- Remarks on surveys where appropriate. 

- An assessment of the adequacy of the survey coverage. 


For those not familiar with hydrographic field sheets, a sample sheet, 
with the derived chart, is included at the end of this introduction. 


Tidal 


The location of the nearest tide gauge, the period over which tide 
heights have been analysed, and the date of analysis. 


Oceanography 


A narrative discussion of the information available. 


Bibliography 


An alphabetical listing by author of the known oceanographic publica- 
tions on the estuary region. 


The hydrographic section of the report was prepared by B.M. Lusk, the tidal 


information was collected by W.J. Rapatz, and the oceanographic bibliography 
was compiled by E.W. Marles. 
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ESTUARY LOCATION ON THE BRITISH COLUMBIA COAST 


Fraser River 

Indian River - Indian Arm 

Squamish River - Howe Sound 

Cowichan River - Cowichan Bay 
Chemainus River 

Nanaimo River 

Puntledge - Courtenay River - Comox Harbour 
Campbell River 

Salmon River - Kelsey Bay 

Homathko River - Bute Inlet 

Quatse River - Port Hardy 

Wannock River - Rivers Inlet 

Bella Coola River - North Bentinck Arm 
Kitimat River 

Skeena River 

Gold River - Muchalat Inlet 


Somass River - Alberni Inlet 
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FRASER 
Hydrography 
1. Surveyed most recently in 1967/68. 


Zea CAL LOL WF pence lols) 


oo 7a) 1:30,000 Strait of Georgia 


F.S. 2272-L 1:15,000 Burrard Inlet 
3. The Fraser River itself is continually being surveyed for depths 
by the Department of Public Works. The last hydrographic survey 
recorded on the river itself was in 1929, sheets 2209-L and 2210-L. 
Scale (12,000. 


4. Chart scale of the river - Charts 3488 - 1:25,000 
3489 - 1:18,000 


5. Chart scale of Str. of Georgia 3450. - 1:80,000 
1:50,000 


Sounding Coverage 


A modern survey with good coverage of tidal estuary. 


Tidal Information 


For a one year period in 1958-59 water level gauges were operated to 
Tidal Survey standards for a 12 month period at the following locations: 


Steveston Woodward Landing Burr Landing 


Fraser Street Bridge New Westminster Port Mann 
(North Arm) 


Pitt River Pitt Lake Hammond 
Wannock Mission Sardis 


A numerical model of the Fraser estuary has been developed, in which 
the Strait of Georgia tides and river runoff can be input and water 
levels at any point in the estuary can be forecast. Several temporary 
tide gauges were operated in the Fraser estuary to calibrate this 
model. 


Oceanography 


The Fraser River Estuary, because of its size, proximity to populated 
centres and importance to the salmon fisheries of British Columbia, is 
one of the few intensively studied estuaries on the Pacific coast. The 
whole Strait of Georgia is often considered to be the Fraser Estuary due 


ma 


to the magnitude of this great river's influence. Oceanography in the 
region immediately seaward of the delta front has been studied by the 
University of British Columbia, Fisheries Research Board of Canada, and 
the Marine Sciences Branch of the Department of Environment over many 
years. A study of the salt water intrusion in the estuary was published 
by Tabata and LeBrasseur (1958). Further studies in the lower river 
have been undertaken in connection with the disposal of sewage from the 
City of Vancouver. Ages (personal communication) and Bawden et al (1972) 
have collected some information on circulation and water characteristics 
on the shallow banks of the delta. 
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INDIAN (BURRARD) 
Hydrography 
1. Surveyed most recently in 1954. 
2. Scale of F.S. 2255=Ses 12245320, 
3. Chart number and scale 3435 - 1:24,300. 
Sounding Coverage 


Very poor bathymetric definition of estuary at this scale. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Deep Cove 1 month 1964 
Lake Buntzen 1 month Y9i2 


Oceanography 


Gilmartin (1962) describes the data record in Indian Arm as "the most 
complete set of oceanographic data available on a British Columbia 
fjord". He has "established in some detail the characteristic annual 
cycle of fjord oceanographic properties". The Indian River influences 
the estuarine circulation within this inlet as a Significant fresh 
water source. A detailed study of the estuary near the Indian River 
mouth has not yet been published but the understanding of the whole 
inlet as an estuarine body of water is much better than average. 
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Outfall flow effect on water quality in the 
Harbour. Rept. no. NV, 750: 30 p. (unpub.) 
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Pickard, G.L. 1961. Oceanographic features of inlets in the British 


Stockner, 


Columbia mainland coast. J. Fish. Res. Bd. Can., 18(6):..907- 
999, 


J.G. - Indian Arm T/S, DO, Zoopl. and P. Prod. Data on 
file. Environment Canada, Pacific Environment Institute. 


Swan Wooster CBA. 1970. Burrard Inlet crossing report K4019/1. 


Tabata, 5. 


Waldichuk, 


Waldichuk, 


T.W. Beak Consultants Ltd., Vancouver, B.C. 


1970. Report on hydrological studies, Stage 2, field obser- 
vations. Nat. Harbours Bd., Burrard Inlet Crossing Rept., 
no. 29, Vol. 1-67 


1970. A brief oceanographic description of the waters of 
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SQUAMISH 


Hydrography 


1. Surveyed most recently in 1939 small scale. 
1930 large scale. 


2. ,e€alesOf-F Set 27s=Sr-Viw 12 ton: 
3. Chart number and scale 3586 - 1:37,500. 


4. WM. J. STEWART will be carrying out a hydrographic survey of 
Squamish at a scale of 1:8,000 in May 1973. 


Sounding Coverage 


Bathymetry in estuary poorly defined at this scale. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Squamish 1 year 1961 


Oceanography 


The recent study of the Squamish Estuary in relation to harbour develop- 
ments, (Dept. of Environment, 1972) provides the first "study in depth" 
of all facets of estuarine ecology and oceanography on the British 
Columbia coast. Squamish therefore joins the select group of regions 

on the Pacific coast with a data base of some substance due to their 
importance in the industrial economy of the Province of British 
Columbia. 
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COWICHAN 
Hydrography 
1. Surveyed most recently in 1960. 
emcee. oF PY6A 1177 -L4-. 112',160% 
3. Chart number and scale 3470 - 1: 18,000. 
Sounding Coverage 


Fairly good bathymetric definition of estuary. 


Tidal Information 


Period’ of Year 700 
Nearest Gauge Analysis Analysis 
Cowichan Bay 1 month 1959 


Oceanography 


Although no studies of the Cowichan Estuary have been published, some 
surface drift information and temperature/salinity data have been 
collected within Cowichan Bay and river mouth. (Herlinveaux, personal 
communication). Surface drift data and temperature/salinity relations 
were also investigated in conjunction with the chinook and coho project 
of the Nanaimo Biological Station in the Cowichan Estuary. Sparrow 


(personal communication) has a manuscript report in progress on these 
data in his files. 


Bibliography 


Herlinveaux, R.H. . Oceanographic features of Cowichan Bay - sur- 
face movements and temperature/salinity distributions. In 
progress. Environment Canada, Marine Sciences Directorate. 


Sparrow, R.A.H.* . Cowichan Bay and estuary, May-July 1967, April 
6-June 14, 1966. T/S profiles and surface currents - with 
H. Godfrey. Data on file. Nanaimo Biol. Sta.* now with Fish 


and Wildlife, Province of British Columbia. 
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CHEMAINUS 
Hydrography 
1. Surveyed most recently in 1962. 
2. Scale of #.S.€1197-S = 1:9,120: 
3. Chart number and scale - 3471 - 1:12,000 
Sounding Coverage 
Fairly recent survey but not too much definition of estuary bathy- 


metry. 


Tidal Information 


Period of Year of 

Nearest Gauge Analysis - Analysis 
Chemainus 1 month 1961 
Crofton 1 year 1961 


Oceanography 


There is again, no published data dealing specifically with the Chemainus 
River Estuary, but Waldichuk et al (1968) have obtained oceanographic 
data in Stuart Channel and along the Chemainus River delta front in 
studying the effects of a pulp mill at Crofton. Fulton et al (1967, 
1968, 1969) have collected some information at a repetitive station in 
mid-channel off the river mouth. 


Bibliography 


Fulton, J.D., 0.D. Kennedy, K. Stephens, and J. Skelding. 1967. Data 
record. Physical, chemical and biological data, Strait of 
Georgia, 1966.) -Fish. Res. Bd. )GaneeMsi Rept. ,2mow91S, 


1968. Data record. Physical, chemical and biological data, 
strait of Georgia, 41967 48 Fish*) Res. wBdy sGangi Ms. Rept), no. 
968. 


Fulton, J.D., 0.D. Kennedy, J. Skelding, and K. Stephens. 1969. 
Physical, chemical and biological data, Strait of Georgia, 
1968. Fishs Res’«Bdy«Can. Ms: Répt:...ino. *1049. 


Pacific Oceanographic Group. 1954. Physical and chemical data record, 
Strait of Georgia, 1949-1953. Fish. Res. Bd. Can. Ms. Rept., 
Nanaimo Biol. Sta. 
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Waldichuk, M. 1955. Effluent disposal from the proposed pulp mill at 


Crofton,+B,C. Fish, Res..’Bd.. Can. Prog akept. Pactiite Coase 
Sta. , no. 102% 6-9. 


1958. Summer oceanography in Osborn Bay, B.C. Fish. Res. Bd. 
Can,” Prog. ReptatPacific Coast. otac,enOw Lh 6-1 on 


1964. Dispersion of Kraft mill effluent from a submarine 
diffuser in Stuart Channel, British Columbia. J. Fish. Res. 
Bd, "Canw, 2h(o) 2) 2269-1516; 


Waldichuk, M., J.H. Meikle, and J.R. Markert. 1968. Physical and 


chemical oceanographic data from the east coast of Vancouver 
Island, 1954-1966. Fish. Res. Bd. Can. Ms. Rept., nos?9s9, 


Waldichuk, M. . Data on file. Stuart Channel 4-17 Aug., 1969 - 


26 July-8 Aug., 1971. Pacific Environment Institute, 
W.) Vancouver, ~B.C. 
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NANAIMO 


Hydrography 


1. Surveyed most recently in 1965/67 - Estuary 
1937 - Harbour 


2. DCOLEROL Or .See 1951, <7 1 522; .000 
1138-L - 1:6,080 


3. Chart number and scale 3581 - 1:6,100 
3456 - 1:18,000 


4. New survey of harbour area scheduled for 1974. 
Sounding Coverage 


_F.S. 1193-L shows low water in the estuary, however no sounding 
was taken during the years 1965/67. 


F.S. 1138-L shows soundings north of this low water line. Sound- 
ings acquired years ago. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Nanaimo 1 month 1926 


Oceanography 


Tully and Waldichuk (1953) produced some surface current information 
below the delta front of the Nanaimo River Estuary. Waldichuk et al 
(1968) obtained some data in Northumberland Channel adjacent to the 

pulp mill established there. The majority of this data is in a restrict- 
ed channel to the south of the Nanaimo River and not directly related 

to the estuary. Some temperature/salinity data is on file at the 

Nanaimo Biological Station. kK. Stephens (personal communication). 

The total volume of data in the Nanaimo Estuary seems to be very small. 


Bibliography 


Henry, R.F., and T.S. Murty. 1971. Three dimensional circulation in 
a stratified bay under variable wind-stress. Third Liege 
Colloquium on Ocean Hydrodynamics, Univ. of Liege, 3-8 May, 
LOL; 


Stephens, K. . T/S, DO: Datavon file, Nanaimo Biol. Sta. 
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Tully, J.P., and M. Waldichuk. 1953. The oceanographic phase of the 
Nanaimo sewage problem. Pacific Oceanogr. Group Ms., Nanaimo 


Biol. sta, 

Waldichuk, M._ . Departure Bay study, 1968. Ms. report in pro- 
gress. Pacific Environment Inst., West Vancouver, British 
Columbia. ; 


Waldichuk, M., and J.P. Tully. 1953. , Pollution, study in Nanaimo 
Harbour. Fish. “Res. /Bda. Can: Prog. Rept. Pacitic Coasts ota, 
noy 972 14-17. 


Waldichuk, M., J.H. Meikle, and J.R. Markert. 1968. Physical and 
chemical oceanographic data from the east coast of Vancouver 
TsYand, 1954-1906. Vol. 2." Fish. Kes. Bd. Can. Ms, kept... 
no. ‘989. 
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PUNTLEDGE - COURTENAY 
Hydrography 
1. Surveyed most recently in 1969. 
2. Scalexso£.F 8.0. 201-L..-.1-.405000 
3° “Chart: numberyand scale"3599 »-),1:'10,, 000: 
Sounding Coverage 
Bathymetry in estuary good and at a large scale. 


This F.S. and Chart used as an example of estuary sounding coverage. 


Tidal Information 


Period of. Year of 
Nearest Gauge Analysis Analysis 
Comox lsyéear 1968 


Oceanography 


Oceanographic data by Waldichuk et al (1968) and Crean and Ages (1968) 
was collected within Comox Harbour. The Pollution Control Branch of 
the Province of British Columbia is operating an investigation of 
oceanographic conditions for water quality evaluation. (Webster, per- 
sonal communication). There has been some work in the Puntledge River 
in connection with dam construction which might bear upon the estuarine 
situation. The oceanographic information does not extend upstream be- 
yond the delta front. 


Bibliography 


Crean, P.B., and A.B. Ages. 1968. Oceanographic records from twelve 
cruises in the Strait of Georgia and Juan de Fuca Strait. 
Dept. Energy, Mines and Resources Ms. Rept., Marine Sciences 
Branch, \Victorilas > vols: 


Waldichuk, M. 1962. Oceanographic characteristics of Comox Harbour 
and approaches in relation to sea disposal of sewage. Fish. 
Res. Bd. Can., Nanaimo Biol. Sta.: 49 p. Mimeo. 


Waldichuk, M., J.H. Meikle, and J.R. Markert. 1968. Physical and 
chemical oceanographic data from the east coast of Vancouver 
Island, 1954-1966. Fish. Res. Bd. Can., Ms. Rept., no. 989. 


Websters... . T/S, DO for water quality. Prov. British Columbia 


Pol l.-Contr.-Br.«e—Data+onv file: 
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Waldichuk, M. . Baynes Sound - Comox Harbour, data on file - 
4-17 Aug., 1969 and 26 July-8 Aug., 1971. Pacific Environment 
Institute, W. Vancouver, B.C. 
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CAMPBELL 
Hydrography 
1. Surveyed most recently in 1950. 
4. pstale of Fas. U16G4=le- gels; 240: 
35. Chart number and scale .3556 - 1:18,247. 
4. Patch on Chart 3556 recently surveyed by revisory surveys 1972. 
Sounding Coverage 


Fairly poor estuary definition at this scale. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Campbell River 1 year 1967 


(permanent gauge) 


Oceanography 


No information directly related to the Campbell River Estuary has come 
to light. Huggett (personal communication) has, on. file, the results 

of some drift pole current studies in Discovery Passage. Waldichuk et 
al (1968) has provided some data for the area north of Campbell River 

from studies in Duncan Bay adjacent to the Elk Falls pulp mill. 


Bibliography 


Huggett, W.S. . Drift pole current measurements in Discovery 
Passage, 1966. Data on file. Environment Canada, Marine 
Sciences Directorate, Victoria. 


Waldichuk, M., J.H. Meikle, and J.R. Markert. 1968. Physical and 
chemical oceanographic data from the east coast of Vancouver 
Island, 1954-1966. Fish. Res. Bd. Can. Ms. Rept.,’no. 989. 
(Several stations off the river mouth.) 


Waldichuk, M. . Alberni Inlet, Muchalat Inlet - Quatsino Sound and 
Discovery Passage Survey, 19 Aug.-1 Sept., 1968. Pacific 
Environment Institute, W. Vancouver, B.C. 
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SALMON (KELSEY BAY) 
Hydrography 
1. Surveyed most recently in 1938. 
Zo esrelesor Ris. 22159. ~ 1:24 ;320. 


3. Chart scale and number 3567 - 1:37,500 
$523 = 14257000 


4. Numerous revisory surveys and patches from engineering drawings 
Since original surveys. 


Sounding Coverage 


Survey at this scale does not show much of the estuary bathymetry. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Kelsey Bay 1 month 1916 


Oceanography 


No oceanographic information in the region of the Salmon River Estuary 
has been discovered. 
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HOMATHKO (BUTE) 
Hydrography 
1. Surveyed most recently in 1957. 
Cteeecdle son FS... 2252-L) = 1:56,480.. 
BP Chart. scaie. ana mumber Oocd =, 1: 75,000, 
Sounding Coverage 


Bathymetry not too well defined at the charted scale. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis Analysis 
Waddington Harbour 1 month 1909 


Oceanography 


Many oceanographic stations have been occupied in Bute Inlet by personnel 
from the Institute of Oceanography, University of British Columbia. The 
repetitively sampled station (no. 8) at the head of Bute Inlet in 
Waddington Harbour is the closest station in this series to the Homathko 
River mouth. This sampling programme was continued from 1951 to the 
present, providing a good time series of observations. Tabata and 
Pickard (1957) have published a study of the oceanography of Bute Inlet 
considering the whole inlet to be the estuary of the Homathko River. 
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QUATSE (PORT HARDY) 


Hydrography 
i 


ra 


De 


4, 


Surveyed most recently in 1942. 
Stale ‘orppes.. 1152-L = 1°12°160, 


Chart number and scale - 3572 - 1:37,500 


insert = 1:12.5200 


Numerous corrections to chart from revisions and public works plans 
Since original survey. 


Sounding Coverage 


Bathymetry not well defined in estuary. 


Tidal Information 


Oceanography 


Period of Year of 
Nearest Gauge Analysis Analysis 
Port Hardy loyear 1970 


(permanent gauge) 


The 1971 float and dye studies (Boyd, personal communication) in the 
Quatse River Estuary appear to be the only work available inside Hardy 


Bay. 


The study was initiated as a result of the proposal by Dist. of 


Port Hardy to establish a sewer outfall at the head of Hardy Bay. 
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1959. Oceanographic data record, Coastal Seaways Project, 
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WANNOCK (RIVERS INLET) 
Hydrography 
1. Surveyed most recently in 1936. 
2. Scale. of P.S/4A22-Lf- 1636, 480. 
3. Chart scale and number 3778 - 1:73,600. 
Sounding Coverage 
Bathymetry in estuary not well defined at this scale. Rather 


ancient. Few revisions since original survey. 


Tidal Information 


Period of Year of 
Nearest Gauge Analysis - Analysis 
Wadhams 1 year 1910 


Oceanography 


Data in Rivers Inlet adjacent to the Wannock River seems extremely 
sparse considering this river's importance as a salmon producing source. 
Fisheries studies in the river and above in Owikeno Lake have been con- 
tinuing for many years. Echo soundings have been collected almost 
daily from 1968-1972 and returns interfering with fish counting have 
been attributed to density structure in the water column (Wood, F.E.A., 
personal communication). Some temperature/salinity information was 
obtained when the echo sounding interference was studied but may have 
been discarded. Aerial photography of plume patterns in the inlet were 
undertaken to investigate gillnet location and effectiveness. 
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Anon 74951.) (British. Golumbia inlet study,’ 195]. Datasgreport..no. “). 
Inst. Oceanogr., Univ. British Columbia. 


1956. British Columbia inlet cruise, 1956. Data report no. 
8. Inst. Oceanogr., Univ. British Columbia. 


Pickard, G.L. 1961. Oceanographic features of inlets in the British 
Columbia mainland coast. J. Fish. Res. Bd. Can., 18(6): 907- 
999. 


Wood ek 2E. A. . Density structure interpreted from echo soundings 
1968-1972 approx. daily. Data on file. Fisheries Service, 
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schutes, 0, . Current patterns from aerial photos of river plume. 
Data on file. Fisheries Service, Environment Canada. 
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BELLA COOLA 
Hydrography 
1. Surveyed most recently in 1956. 
2... Scalervor fo. 4495-8 - 1318 7240. 


3. Chart number and scale 3730 - 1:75,000 
‘Whart =" 1229500 


Sounding Coverage 


Bathymetry good at a good scale. 


) pete ait orma eso 


Period of Year of 
| Nearest Gauge Analysis Analysis 
) Bella Coola 1 year 1966 


Oceanography 


Oceanographic information in the Bella Coola Estuary does not extend 
above the delta front. Dodimead and Herlinveaux (1968) discuss the 
region immediately below the river but the bulk of their data is 
| further down the channel, in North Bentinck Arm and in other reaches 
of this inlet. Goyette (personal communication) reports data on file 
including temperature/salinity and dissolved oxygen from one survey. 
The data available cannot adequately portray the variability of oceano- 
graphic parameters in the Bella Coola River mouth. 
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Anon. 1951. British Columbia inlet study, 1951. Data report no. l. 
Inst. Oceanogr., Univ. British Columbia. 


1956. British Columbia inlet cruise, 1956. Data report 
no. 8. Inst. Oceanogr., Univ. British Columbia. 


Dodimead, A.J., and R.H. Herlinveaux. 1968. Some oceanographic 
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Goyette, D.E. . One survey T/S and DO. Data on file. Environ- 
mental Protection Service. 


Herlinveaux, R.H. 1968. Drift card releases and recoveries in Burke 
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Herlinveaux, R.H. . Some oceanographic features of a British 
Columbia inlet, Burke Channel, 1966-1967. In progress. 
Environment Canada, Marine Sciences Directorate. 


Surface water movements in a central British Columbia 
inlet. In progress. Environment Canada, Marine Sciences 
Directorate. 


Pickard, G.L. 1961. Oceanographic features of inlets in the British 
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NOTE Aluminum Co. of 
Aids to Navigation at Kitimat 
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KITIMAT 
Hydrography 
1. Surveyed most recently in 1952. 
2, Scaleyvof Ps S.cG523<Sn= UtT2guG0: 
3. Chart scale and number 3736 - 1:12,165. 


4. Corrections to Chart 3736, as far as can be determined have all 
come from various plans and dredge surveys. 


5. No revisions or further hydrography has been gathered since 
original survey. 


Sounding Coverage 
Bathymetry in estuary is good though construction additions to 


mills etc. have altered original western quarter. 


Tidal Information 


Period of Wear of 
Nearest Gauge Analysis Analysis 
Kitimat 1 year 1953 


Oceanography 


Like many others, the Kitimat River Estuary remains largely unstudied. 
Waldichuk et al (1968) have collected four series of oceanographic 
survey data in connection with the pulp mill established at the head 

of Kitimat Arm. The Environmental Protection Service has undertaken 

two surveys of temperature/salinity and dissolved oxygen, also connected 
with the pulp mill. The University of British Columbia has only one 
data record with information in this inlet. 


Bibliography 


Anon. 1951. British Columbia inlet study, 1951. Data report no. 1. 
inst Oceanogr., Univ, “British Columbia. 


Goyette, Der. . Two surveys of T/S and DO in Kitimat Arm. 
Environmental Protection Service. Data on file. 


Pickard, G.L. 1961. Oceanographic features of inlets in the British 
Columbia mainland coast. J. Fish. Res. Bd. Can., 18(6): 907- 
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Waldichuk, M., J.R. Markert, and J.H. Meikle. 1968. Physical and 
chemical oceanographic data from the west coast of Vancouver 
Island and the northern British Columbia coast, 1957-1967. 
Fish. Res. Bd. Can. Ms. Rept., no. 990. 


Waldichuk, M. - Fisher Channel - Cousins Inlet, Douglas Channel - 
Kitimat Arm and Prince Rupert area oceanographic surveys, 
15-28 Sept., 1969. Pacific Environment Institute, W. 


Vancouver, B.C. 


Ocean Falls, Kitimat and Prince Rupert, 17-30 July, 
1972. Pacific Environment Institute, W. Vancouver, B.C. 
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SKEENA 
Hydrography 
1. Surveyed most recently in 1949. 


2s * ocalée Givr.S, 33 18=L 
S51/=L 


1:18,240 ) 


1:48 2400 


3301-L - 12245320.) 
3302-1, — 1:24 320 jue 


3. Chart number and scale 3756 - 1:23,100 
DTS SF 1231.700 


4, WM. J. STEWART will be conducting a hydrographic survey in the 
Kennedy Island area at a scale of 1:30,000 in 1973. 


Sounding Coverage 


3318-L and 3317-L - Good definition of river and estuary but now 
24 years old. 


3301-L ) 


3302-L ) Good definition of river but now 59 years old. 


Tidal Information 


Tide gauge stations were operated for several months in 1949 at Khyex 
Point and Kwinitsa. Gauge stations were previously operated at 
Claxton, Port Essington, and in the Ecstall River. 


Oceanography 


There is little oceanographic information in the shallow inshore region 
of the Skeena River Estuary. The influence of Skeena River water ex- 
tends some distance offshore. The work of Cameron (1948) and Trites 
(1956) is concerned with fresh water entering Chatham Sound. Waldichuk 
et al (1968) studied the Wainwright Basin - Porpoise Harbour area in 
relation to pulp mill pollution and took several stations in the Skeena 
River mouth. The bulk of this information refers to the restricted 
passages north of the estuary. Higgins and Schouwenburg (1973), in a 
continuing programme, have obtained some temperature/salinity data in 
the estuary in conjunction with Fisheries studies. The Hecate project 
is a source of considerable data offshore from the Skeena River. 
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Cameron, W.M. 1948. Fresh water in Chatham Sound. Fish. Res. Bd. 
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Herlinveaux, R.H. . Time series T/S data on file by W.M. Cameron. 


1948. Environment Canada, Marine Sciences Directorate. 


Higgins, R.J., and W.J. Schouwenburg. 1973. Biological assessment of 


Mackay, B. 


fish utilization of the Skeena River Estuary with specific 
reference to port development in Prince Rupert. Fish. and 


Marine Service Tech. Rept., 1973-1. 


S, 1954. Tidal current’ observations in Hecate Strait, J. 
Fish. Res. Bd. Can., 11 (2 )se48-56:. 


Pacific Oceanographic Group. 1955a. Physical and chemical data 


record, Hecate Project, 1954. Nanaimo Biol. Sta. Ms. 


1955b. Data record. Current measurements, Hecate Project, 
1954. Nanaimo Biol. Sta. Ms. 


1955c. Physical and chemical data record, Hecate Project, 
with appendix 1, current observations, 1955. Nanaimo Biol. 
Sta. Ms. 


1956. Physical and chemical data record, Dixon Entrance, 
Hecate Strait and Queen Charlotte Sound, 1934, 1937, 1938, 
1951. Nanaimo Biol. Sta. Ms. 


Trites, R.W. 1956. The oceanography of Chatham Sound, British 
Columbia. J.°Fish; Res 77Bd. Can. , 13(3)2685-454., 

Waldichuk, M., J.R. Markert, and J.H. Meikle. 1968. Physical and 
chemical oceanographic data from the west coast of Vancouver 
Island and the northern British Columbia coast. Fish. Res. 
Bd. Can. Ms. Rept... ,. no...990.. 

Waldichuk, M. . Unpublished data record. Oceanographic studies 


in Prince Rupert area, July-August, 1967. Pacific Environment 
Institute, W. Vancouver, B.C. 


Fisher Channel - Cousins Inlet, Douglas Channel - 
Kitimat Arm and Prince Rupert area surveys, 15-28 September, 
T9O9 ) oro! WN aeVanGouver, —B oC. 


Ocean Falls, Kitimat and Prince Rupert, 17-30 July, 
1972. PEL, West Vancotiver .B..C. 
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GOLD (VANCOUVER ISLAND) 
Hydrography 
1. Surveyed most recently in 1967. . 


2. Scale of F.S. 1240-S -.1:10,000 
1200-L - 1:1,200 Wharf plans 


3. Chart number and scale 3665 - 1:12,000. 
Sounding Coverage 
No estuary bathymetry on F.S. 1200-L. 


Very little bathymetry of estuary on 1240-S. 


Tidal: Information 


Period of | Year of 
Nearest Gauge Analysis Analysis 
Gold River 1 month 1967 


Oceanography 


Although there is no published data dealing specifically with the 
estuary of the Gold River, the sampling locations of Waldichuk et al 
(1968) and those of the Institute of Oceanography, University of British 
Columbia, are located in the main channel of Muchalat Inlet as shown 

on the accompanying chart. Oceanographic information in this area was 
collected to evaluate conditions for effluent disposal from a pulp mill 
established at Gold River. Several series of temperature/salinity data 
have been collected by Environmental Protection Service personnel and 
are on file. Goyette (personal communication). 
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SOMASS (ALBERNI) 


Hydrography 


1. Surveyed most recently in 1959 - 1172-L. 
after Tsunami 1964 - 1225-S. 


2. socale of F.S,. Li72-L, = >1s12,760 
1225-5 = P6080 


3. Chart number and scale 3609 - 1:38,800 
Insert - 1:9,000 


Sounding Coverage 


Good bathymetric definition of estuary. Numerous revisions to 
area since original surveys. 


Tidal, iInrormation 


Rer vod! of Year of 
Nearest Gauge Analysis Analysis 
Port Alberni 1 year 19.78 


(permanent gauge) 


Oceanography 


The Somass River Estuary is one of the more heavily studied areas of 

the British Columbia coast because of a pulp mill established in Port 
Alberni discharging effluent at the mouth of the Somass River. There 

is extensive oceanographic information in the area directly below the 
delta front of the river. A baseline inventory of physical and bio- 
logical conditions of the delta region has been prepared. (Howard 

Paish and Assoc. Ltd., 1972.) The data records in the bibliography 
include cross seasonal and annual coverage for long terms and therefore 
could provide one of the best data bases available in a British Columbia 
estuary. 
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1960. Waste disposal in the marine environment. Pergammon 
Press 


Waldichuk, M. 1956. Pulp mill pollution in Alberni Harbour, British 
Columbia. Sewage and industrial wastes, 28(2): 199-205. 
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INTRODUCTION 


In 1971 the Canadian Hydrographic Service development group, under 
A.R. Mortimer, tested the Trisponder Model 202A pulsed radar positioning 
system. Some changes have been made to the system since then, the most 


significant being new transponders and antennae. 


In September 1973 Computing Devices of Canada Limited lent the system with 
new Model 210 transponders to the Hydrographic development group for range/ 
accuracy evaluation. During the few days the system was available signal 
stability tests were conducted on a survey launch over a distance of 30 kilo- 
metres (km), and precision and stability tests were conducted over seven 

land baselines previously established by tellurometer traverse and triangu- 


lation (Anderson et al, 1973). 


Following these tests, the system was evaluated by the B.C. Forest Service 
for positioning a helicopter used for aerial photography. Mr. Woods assisted 
in the airborne operations and further observed the system's operational 


characteristics. 


THE EQUIPMENT 


A. The Distance Measuring Unit (DMU) 


This unit is basically the same as that tested by A.R. Mortimer (1972). 
The DMU contains the range measuring circuitry and operating controls. 
Two digital displays present the measured distances in metres (m). Its 


power consumption is 40 watts and it requires a 22-32 volt d.c. power 


supply. 


The system measures the distance to each remote transponder by using an 
X-band frequency and averaging 10 or 100 (depending upon front panel 
control setting) sequential measurements. An average of 10 readings 
resolves the distance to the nearest 10 m and 100 to the nearest metre. 
Remote identification is achieved by utilizing separate coded pulse 
repetition frequencies (PRF) for each transmitter. The averaged meaSdure- 


ment is displayed on the DMU digital display at one second intervals. 


The DMU has 4 separate calibration screws for 4 transponders, any 2 of 


which may be selected for use at any one time. 


The model tested was equipped with a ''time sharing option" that enables 
up to 4 master units to operate simultaneously with the same 2 trans- 
ponders. It also features a manual update control, in which the last 
range measurement is held until the operator manually initiates a new 
reading. This feature may be particularly useful during shoal exami- 


nations. 


Red fault lamps warn of signal loss or interruption during the measure- 
ment sequence which may be caused by a loss of line of sight, loss of 


transponder power, etc. 


The Master Transmitter/Receiver (Master Transponder) 


The master transmitter interrogates the remote transponders and receives 
the transponder's reply. It utilizes a 6 decibel (db) omnidirectional 
antenna which is coupled to it by a simple quick release connector. The 


; : é O 
antenna's vertical beam width is 30 . 


The transmitter's magnetron is tunable from 9300 to 9475 megahertz (MHz), 
and its power output is 1000 watts at +60 dbm (power output in db above 
1 milliwatt). Power consumption is 17 watts. The receiver sensitivity 
is -74 dbm (sensitivity in db below 1 milliwatt). The unit requires a 


22-32 volts d.c. power supply, which is supplied by the DMU. 


The Transponders 


The transceivers of the remote transponders are identical to the master 
transceiver. However, instead of an omnidirectional antenna, they use a 
16 db vertical slotted array directional antenna with a horizontal beam 


width of 87° and a vertical of te 


The transponders operate from 22-32 volts d.c., which is normally sup- 
plied by two 12 volt batteries, and the power consumption is 17 watts. 
Forty minutes after interrogation ceases, the transponders switch to a 


standby mode in which the power consumption is reduced to 8 watts, thus 


a 


extending the period in which the remote stations can be left unattended. 


All transponders can be mounted on standard tripods, European tripods, 
or a one inch threaded pipe. The manufacturer's specifications are 


listed in Appendix I. Plate 1 illustrates the equipment. 


INSTALLATION AND OPERATION 


The portability of the DMU allows it to be installed in any convenient place 
aboard a launch, although it is not adaptable to a standard 19 inch rack. 
Fifty feet of cable is supplied to connect the master transceiver to the DMU; 
since the antenna and master transceiver are one unit, this allows the anten- 
na to be attached to a mast or other suitable location. During the launch 
tests the master unit was mounted on a Wild tripod lashed to the mast step. 
The remote transponders, also mounted on tripods, were located at known shore 


stations. 


Before operation, the system must be calibrated over a known distance; pref- 
erably this distance should be one approximating the central distances used 


during the survey. 


Since the DMU provides binary coded decimal (BCD) output, peripheral equip- 
ment may be used aboard the sounding vessel, such as a remote read-out, track 


plotter, left/right indicator, etc. 


METHOD OF TESTS 


A; “Calihration 


At the calibration distance of 2102.4 m, transponder A was reading 1 m 
high and B was 1 m low. Rather than adjusting the read-out on the DMU, 


the appropriate correction was made to the raw data. 


B, ‘Stability Aboard a Survey Launch 


In order to determine the stability of any one distance reading, the same 


tests were conducted as had been used for the Mini-Ranger tests in 1973 


a4 


and the Trisponder tests in 1971. The method is to put both transponders 
at the same shore station, approximately six feet apart on a line 90° to 
the launch so the distance from the launch to each transponder is essen- 
tially the same. The launch is then steered along an arc of constant 
distance at the centre of the beam pattern and both range readings are 
recorded simultaneously. The readings are then plotted on a graph and a 
curve of best fit is drawn to eliminate erratic movements of the launch. 
The dispersion of the readings about this line, plotted on a straight 


line graph, show the measure of instability, (see Figure 1). 


In our tests launch breakdown and lack of time allowed for just a single 
run at 30 km. The root mean square (r.m.s.) of the dispersion of read- 
ings about the actual launch course is +1.7 m, Stability at shorter 
ranges should be the same or better unless the system is affected by out- 


side interference. 


Stability of the Trisponder may be affected by "ducting", similar to 
other pulsed radar positioning systems (Anderson et al, 1973). While 
steaming straight away from the transponders the readings began fluctu- 
ating at 15,500 m and continued until approximately 16,000 m. Stability 
broke down even to the point where the radiofrequency (RF) link to the 
transponders was broken and the "remote off" fault lights flickered on 
and off. It was at about this same distance at the same location, that 
the Mini-Ranger experienced stability breakdown during the tests in the 
spring of 1973. Such occurrences are random in nature and may or may not 


happen in various geographic locations. 


Stability of the System Over Precise Land Baselines 


Land tests were conducted with the master unit set up on the Mount 
Douglas Geodetic station and the transponders were taken by car to seven 
previously measured points from 2 to 21.5 km, These are shown in Figure 
2. At each station the transponder antennae were pointed directly at the 
master (0°, or the centre of the beam pattern) and 10 readings were 
recorded. The transponders were then turned to 20°, 40°, 60°, 80°, 90°, 
100°, 120°, 140°, 160°, and 180° from the centre of the beam pattern 


enabling the accuracy to be examined at the edge of the beam pattern. 


Bo. 


Unfortunately, some data for transponder B had to be deleted due to 


pointing errors caused by fog. 


The dispersion (or standard deviation) of the 10 readings taken at each 
angular increment may also be regarded as the measure of stability. The 
most instability was found at a distance of 15.4 km with the other ranges 
showing stability equal to or better than this range. Therefore, by com- 
bining the values at 15.4 km and 21.5 km a value for the entire test 
range is found. The root mean square (r.m.s.) value at 15.4 km, centre 
of the beam pattern, is +1.8 m and at o.. it Sis te te, At 21.5 km the 
r.MpS, £07, Stability, at_0° is.+1.2 m,and_at 90° it is +2.1 m. By combin- 
ing these four figures a total r.m.s. value for. the stability of the 


system from 2 to 21.5 km and covering 90° of the beam Dattermnaasy +Zad. mi 


In order to determine instability caused by roll or pitch of a moving 
launch; the master antenna was tilted towards and away from the remote 
station (which was pointed directly at the master) in increments of 10°. 
Signal instability occurred at 30° from the vertical which should not 


cause any problems under normal sounding circumstances. 


D. Precision 


The. precision of any one position line within the beam pattern of the 
transponders can be determined from the mean differences of the range 
readings from the known baselines. These differences from each known 
distance are tabulated in Appendix II. The total r.m.s. value about the 
median of the mean differences, from 2 to 21.5 km is +6.5 m, between 

0° and 90°. Figures 3 and 4 show a graphical representation of the mean 
differences for the system calibrated at 2 km. Figure 5 shows the pre- 
cision contour for the average of all readings. As with other radar 
positioning systems, these contours can be changed by altering the dis- 


tance at which the transponders are calibrated. 
TOTAL SYSTEM ACCURACY 


In the previous Trisponder evaluation stability and precision were combined 


with a third factor, repeatability, to give total system accuracy. Since we 


say 8 


had the Model 210 for only a few days a series of calibration checks could 
not be made to determine frequency drift. However, if we assume that fre- 
quent calibrations made over a long period of time would maintain steady 
repeatability, we can use +1 m. This assumption is based on the original 
findings of A.R. Mortimer (1972) and the fact that the DMU is essentially 


the same unit. 


Therefore total system accuracy from the land tests is: 


2 2 2 
Total = Ey + E, + ER 


leo eal + eons? + +12 


+6,.9 m (throughout 90° of the beam pattern) 


= precision error, and E, = repeatability error. 


E, = stability error, E 3 


1 2 


Figure 6 shows the overall accuracy contours for the system calibrated at 
2 knis 


RADAR INTERFERENCE 


During the launch tests (conducted in close proximity to several vessels 
operating radar), no radar interference was experienced; a noticeable differ- 
ence from 1971. This is due, apparently, to changes made in the gating tech- 


niques used in signal reception. 


Each transponder receives a series of coded pulses at a particular pulse 
repetition frequency (PRF), each PRF being different for each transponder. 
Signals are sent back to the master at the same PRF as received, thus distin- 
guishing one transponder from the other. Each transponder has an electronic 
"gate'’ which is only open to accept its particular PRF for 2 1/2 microseconds 
out of a 1 millisecond interval (or about .02% of the time), thus eliminating 


outside interference such as radar (von der Heide, 1973). 


HELICOPTER OPERATION WITH TRISPONDER 


Following our tests the Trisponder was lent to the B.C. Forest Service. 

Their application of the system was to position a helicopter during low level 
aerial photography of various stands of timber. This exercise took place in 
the Cowichan Valley near Duncan, B.C. The transponders were placed so that 
line of sight would be maintained in the working area. To check the system's 
Stability a set of readings from one transponder was recorded while the 
helicopter was flown along an arc of 10,500 m. Normally a curve of best fit 
could be applied to this graph but it was so smooth that little or no dis- 


persion due to instability can be detected. This graph is shown in Figure 7. 


Overall, the precision and stability of Trisponder aboard a helicopter should 
be as good as on a launch provided the aircraft stays within the vertical and 


horizontal beam patterns of the antennae. 


CONCLUSIONS 


The full range and reliability of the Trisponder with new "Model 210'' trans- 
ponders has not been evaluated since the system was available for only a few 
days; however, the new transponders seem to have solved the main problems 


encountered with the earlier system evaluated by Mortimer (1972). 


Accuracies deal with one position line only and accuracy breakdown due to 
geometrical configuration will have to be accounted for when using the system 
in a range/range mode. Accuracies of +6 m should be attainable up to a range 
of 30 km as long as the system is calibrated within the working area; for 
example, to survey between 15 and 30 km the system should be calibrated 
between 20 and 25 km. For a range of 12 km the transponders could be set up 
with the centre of the beam pattern perpendicular to the shore and an accu- 
racy of +5 m could be maintained for a distance of 7.5 km along the shore, 


for a calibration within 2-12 km. 
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Figure 3 
Land Based Precision Tests 
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Figure 4 
Land Based Precision Tests 
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Total System Accuracy 
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Trisponder Model 202A-MS02-M2 


The ComDev Trisponder Model 202A is 
an automatic distance measuring system 
which provides a mobile vehicle with a 
highly accurate position relative to two 
fixed locations, up to a range of 15 miles 
line-of-sight (up to 50 miles line-of-sight 
with the long range system). Using radar 
principles and digital techniques for 
measuring the elapsed times of rf sig- 
nals between the mobile master and two 
remote slave stations, the Trisponder 
System accurately determines and simul- 
taneously displays the resultant two 
ranges. The mobile’s position can then 
be trilaterated manually or by automated 
methods depending on which of the op- 
tional peripherals are interfaced with the 
basic Trisponder System. 


System Features 


Minimum size and weight — ideal for small marine In x1 mode the average of 10 (optionally 100) 
craft and airborne applications independent range measurements are displayed 
and updated at one-second intervals 


Solid state electronics built to withstand tropical 


to Canadian arctic environmental conditions Simple range calibration adjustment 
Built-in capability for selecting any two of four Automatic call-up — remote station activated by 
slave transponders interrogation signal; reverts to standby when 


interrogation terminated * 


Low voltage monitor provides automatic shutdown 
before battery failure « Mobile antenna — omni/30° slotted waveguide 
Remote antenna — 87/5° vertical slotted array 


Standard BCD logic outputs 


Master and remote units interchangeable — 
minimum spares required * * Available with long range transponders 


Init Specifications | 
TRISPONDER MODEL 319X 
the short range system 


RF Power 
Frequency 
Environmental! 
Power Supply 
Size 

Weight 
Housing 


luding antenna) 
herproof drawn aluminum case 


the long range system 


Transmitter 

Pulse Power 
Receiver 

Tuning Range 
Antenna Cannector 
Environmental 
Power Supply 
Power Consumption 


Automatic Call-up 
~ Operate 
~ Standby 


Automatic Power Control 
Size (less antenna) 
Weight 


Range 


ange system) 


Accuracy 
Resolution 


Display 
Calibration 
Outputs 


Power Supply 
Size 
Weight 


Housing erproot al uminum transit case 


Options 


100-SUM OPTION Stable readout obtained by summing and averaging 
100 distance measurements from each remote station. 
Display is updated at one-second intervals. 


TIME- SHARING ADAPTOR Allows two or more mobile stations to operate with 
the same pair of remote stations. 


AUTOPLOT SYSTEM Provides conversion of range data to map coordinates 
in Transverse Mercator, Lambert, or other grid 
system for track plotting directly on an X.Y. Recorder. 


ANTENNAS - REMOTE 45°/45° horn for airborne operations. 
87°/5° vertical slotted array for marine operations. 


AUTOMATIC POSITION RECORDING Provides printed record of Trisponder range data at 

SYSTEM APRS- 301 selected intervals of one second to one hour. 
Optional printouts of time and BCD data, e.g., depth, 
temperature, etc., can be included. 


RF LINK SIMULATOR The operation of transponders can be simulated for 
testing the Distance Measuring Unit. RF outputs 
representing a variety of ranges can be individually 
selected. 


PORTABLE RF DETECTOR Received signals are converted to an audio output to 
confirm that mobile/remote stations are radiating. 


LEFT” RIGHT TRACK INDICATOR Enables a mobile station to steer a course parallel to 
any preselected line of position. 


REMOTE DISPLAY UNIT Provides a remote display of range data from the 
Trisponder DMU. 


PN 73022 MAY 1973 PRINTED IN CANADA 


model 210 
microwave transponder 


SReriiieieennnnoiancaionniiiccnnnconaconanansannaninnnnananavoaroonaanaaninannannnncion naanngouneouaaaronnoaaes 


the only. 


Seiiairiianeaiionicananninioiiiooonosnas anne snnnnnovo cio nitanancnatoncsanmanannaian nasa: 


' 


* efficient: requires only 17 watts operating, 8 watts in 
stand-by, with automatic transfer between stand-by 
and operate modes. 


+ 


* accurate: advanced video processor and receiver provide 
= 3 meter accuracy to 80 Km. 


* powerful: a 1 Kw transmitter and —74dbm receiver provide a 
50-mile range with an omnidirectional antenna. 


* reliable: designed for continuous operation; internally 
protected against low voltage, reversed polarity, 
and adjacent high-powered radar. 


portable: engineered and packaged for rugged marine, land, 
and air craft environments— yet easily hand carried. 


* accessible: all components are easily removed. 


Developed for the , Del Norte’s Model 210 is the only transponder designed spe- 
cifically for positioning systems and thus offers significant advantages over general purpose units. 


Ranges in excess of 50 miles are obtained with standard 202A vertical slotted array antennas; ro- 
tating antennas are not required. An advanced receiver and video processor result in nearly con- 
stant pulse delays at signal levels as low as —74dbm. Thus, an accuracy of +3 meters is main- 
tained to maximum range. Performance is further enhanced by a digital range gate which rejects 
interfering signals. 


A highly efficient power converter and extensive use of C-MOS circuitry result in an operating power 
input of only 17 watts. Power used by remote stations during stand-by periods is reduced to 8 watts 
by the Automatic Call-up feature. Call-up from stand-by to operate is initiated by a coded interroga- 
tion from the master station; transfer to stand-by is automatic 40 minutes after normal interroga- 
tion ceases. Based on an 8-hour operating day, a Model 210 remote station will operate 12 days on 
two 12-volt 70 amp-hour storage batteries. 


In designing the Model 210, components with reliable performance records were selected. Circuit 
parameters were so set that operation is well within established ratings. For example, the transmit- 
ter employs a hi-re! magnetron operating at 14 its normal rating. This design approach has resulted 


in a transponder which provides dependable performance in continuous duty service. 


Access to the interior of the transponder is accomplished by removing four self-sealing screws and 
slipping off the waterproof cover. A majority of the electronics are on plug-in PC boards and, as 
potted modules are not used, components may be checked and easily replaced. Since the master 
station is identical to the remote, only one spare is required for complete back-up. An optional test 
fixture allows system gain checks plus transmitter and receiver retuning to be accomplished easily 


in the field. 


The Model 210 is compatible without change in any system; engineering assis- 
tance is available for other applications. 


RGN GO fn 3 iets Wednieeta seh pee 300 feet to 50 miles (100m to 80 Km) 

Transimitferos: snes 5.427. o ae Magnetron, tunable, 9300-9475 MHz 

Pulse Power ................ 1000 watts (+ 60dbm) 

Pulse: Width...) oo nsen e228 0.5 +0.1 microseconds 

Receiver ise oe ns —74dbm, solid state, tunable, 9300-9500 MHz 

Intermediate Frequency ...... 60 +5 MHz with 12 MHz bandwidth 

Power Input: ce ocd Stags Be 22-32 vdc, 8w stand-by, 17w operate 

Automatic Call-up ........... 1 min. by coded interrogation from master 
Stand-by wl 5 fn 4 40 min. after last interrogation 

Size, less antenna ........... TAS hex: WG 1085s di(35:650 I5:ci 2 627 cm) 

Weight-).)-to ware Lies ee 15# (6.8 Kg) 

MOOnTING 56 sists a ce se tmbesor it. U. S. and European Tripod plus 1” NPT 

Antenna Connector .......... UG40/U wave guide flange, quick disconnect 

Wave Guide ............... WG-16, WR-90, RG-52 

Power Connector ............ PTOOE-14-12P, mating connector PTO6W-14-12S 

Temperature ............... —30°C to +-70°C (—22°F to + 158°F) 

Altitude eon oe okie. . 760 to Imm Hg (with vent) 

POCK mC Gy savers oiciee nae > uss Marine environment 


del norte technology, inc. 


P.O. BOX 696 EULESS, TEXAS 76039 (817) 267-3541 


APPENDIX II 


TABULATION OF LAND BASED TESTS 


eat 
— 7 ae Pe 


es oy aes ine as oe 
‘stabenan) tie Nis p- = a 


an jh ainbhy Pi er ere ins! ging © ioe is Jas Se C4 veulity ; 


by sad vai 


wnd.V ahi oo ' ‘y erwape swt wey anit pe ely 
i: hy Oy ears ; og Fie y tov ove seer 


Mice uae : 
"aaali—s e 
0 me ih hee 


a ; 
ins 4 an 
rn: - 


oun 


a nm (here : a ; 1 ely a ‘aygte @h ane wil oP ore (2 a 
a sah. WARE \ pes air | - il ; 
: al | 
; f . byl we TOe CA gnad ware of teenie, 
any ‘i ' j he teryshs Paes Quali File evenviisie, hid tue 
. é ; : : 0 Ged ci.0r 46h Ge 
’ ny Ey) i 


dhe fr a ‘ ee! seating etree 
sae. ’ shad Yeray Pw) he Te ards due, 5 
_, . | neice qm Te PaTPAL ci nee the Moaer 
” : qiiraa bo plete lear k. vy. An optionatitg 
hed ac wy ¥ 


CECE SHEESH CHESKERE FHEEES BL S°de€ 8°2 9°e al | 8°Te WH IZ 
ERSESE KHESHE REKESEH BEEERE YE S°t2 9°T veil €°2 9°¥1t Wy SGI 
o°T t°9 ral | 2°e £0 1°S Sade | Z2°e€ 9°0 ty wn 6 
e°2 rAd tee 2°21 g°0 e°g 2°0 2°*9 5 ode | 6°S wm 2 

& e e e s e 9°2 Ww» 9 
9°0 oe v° tl 2°€ e°O C°zZ 8°*O e°t eae | 2°! Wy 8 
e*o 9°6 2°2 G°St Z*O0 = £=6°*9 ¥°0 z° S$°0 S°e ww Zz 
°a°S NVAW °0°S NV3SK °“aG°S = NWaW °a°S NWSW £°4a°Ss NW3WN 


G°*2.... £°2 S°t 6*ea! Ger. 9°et... at. Ne2t eet 6 °¢€) £°0 OF tt WM Te 
Z°z S*2! t°z2 6° tt wt S°s 1*2 6°9 2°t 0°98 g°2 e°2 Wy GT 
2°*0 1* 6°0 1*2 2*0 S*l 9°0 ret 2°0 2°O 2°0 9°O Ws 6 
£°t °S | ode | c°v 2°y g°2 e°t Set €°t 1*2 zl 7 mn | WM 2 
Ee | ee ee | cen 9°... Gro... ore #*O- 6°0 *D- °1 v*t- 4 9 
e°t €°t e@°T 9°0 Te? 2°O 2°} Te t— o°l 9°T— e°T E°tI- WH 
O° 9°2 t*T G*¢ 9°0 8°O z2°0 S$°O0— ort 1°0 T°T 2°0 wm 2 
eaes Nvaw § °a°s NV SW °0°S NWSW °-a°S NW3W j$“G®S NV3WN £°G°S NVaw 
°s930 02 *s53q0 0 3ONVY 


(S0ueysTp JejyowoINT Te} snutw sSutpees aepucodstay OT jo ese Ieny) 


*O°S = NOILVIAAGD OYVONWIS _. NVWAW = 3ONSe43Sa5dI0 NVA v YSONDdSNVEL 


Pa Ss 


: | aetemow i fad emai aqnibess webnogalsT of to epstevA) 
= eteiG & veesa os se030 6 


waa oQe% WAI +Ge2 WARM oGee WARM 


Rat Evo tet ‘ux s 
—7e S$ --B-- oie 
E«i «S$ Sef Seal “xf 
T.0° 6 -. °Te0 3-0 me 2 
Sei 9.8 Ess Sex hs a : 
S.£ €.6)... ¥ «0 33 
“Suse eT esse 681 as oT a 


SREKRH SEAS KCHEREEA SKEEEE KCHEKKKS SHERRH FZ OF Le. Set €°02 WH T2 
TEKKEE SHEKRK BKEKKSE SHSSKE RKEKKE HSSEKEK HOH GS*°9E | hr 4 O*fc WH ST 
9° 9°sS Se°t ¢°s 2et tT°S ¥*t 2°sS en | €°v Wm 6 
6°0 0°01 y am | S°eT Pde | ¥°9 S°o 2°Z Ade | 1°9 WH 2 

bd ° vet Z2°*9 (eae | 2°sS wa 9 
oO°T 6°E Tad | 1°S wet 6°2 sade | = ae | S°tT i asl WH & 
te? G°SI o°t 2°*ct e*O - w*2Z 9° he & | 2°0 9°S ws 2 
*a°Ss NVoW "aes NV3W *a°S NV 3W ears NVW2W eaes NW 3W 


a ¥°02 £1 S°et Shoes Eire 9°eT. 8°00. e°st 4°90. Lt eee ee. 
e°? 2°! T°? ¥°6 B°t 2°56 wy ST 
B°O €°sS eed | E°e | ie | S°e€ s°! =e eer 0°2 Bet 0°2 Wx 6 
2°0 veg 2°0 v°sS S°0 9° 2°0 GPE ae | B°E 2°0 2° WH 2 
= ee °Oo 6° O° 2°2e ws 9 
6°0 O°2d 6°0 e°T o°T vel | Hd | es o°T v*0 6°0 0°0 WH Yo 
| ae | O° iat | hgh ® ¥°0 Lae | 1St £°0 2°0 T°O0=- 9°0 c°0- WH 2 
*O°*sS NV AW aS SS NVAW 7:0 >>. NVW3W “aes NW SW [a°s NWSw wl fall = NVAW 


*s$930 OF *$s930 02 


(S0ueysTp JeqowoInTTe4 snutW sSutpeet repuodstiy OT jo edeseay) 
"O°S = NOJAVIASO OMVONVAS ONY SW = SONSYaSaSIO NVSW =) YSIGNOGSNVUL 


“a5 


re Sete S28 eet Bei ted i Fee S2€é 4 ° « 
#s $22 Pei -sF a Bas. a : PI 
a eas "Ee all tei Fé Pei eth $s rt: os 
é <I Bat? Sei sok &«9 tet Sei 5<2 wz t 

ae ; ae Toi La? hel S<-f ee Cad a2 = 
Seeunes SGOSGH2 S$24642F £46949 FESesta BHaH EC + od eet ; ded 92S Mp Bi 
= S¢G¢eso8  S9ESFS2 42ESSE 46266 A2SSGUS sé anit ef 2,05 Ws 25 


— a 


eis ; os 7 i pS 


1 a in e, vom Be ; 
aa od ie ney ne i a 


are 
. me: | at 

s ty 
i 


Teepe oe ae NP Ye Tl diate! 


one 
‘eee 


oa eeu Tm 


et ee dais 
PP lt tole = 3 
nee Pn iene aye ee. 
Fare ee ie ie 


on EE 
[i PP IOI TT a sec 


; ce ial menace 

ee nee ey : 
a oeures eens 

Spee ener 3 

: ee a ae papell 


pee 
pan ET ~ 
on < om Sere piace o 
ee ds ¥ E shins ew eer 
See : s sana 
F. : ape yn ae ee 
fs PPE ela IE Pe 


See ni 
a 


: ee NS 
se ae aT eae I LS ST 

ee ence tin ee 
nage i. 


reaper 
Skah aeete r fae 
Teer ee: eens 2 
: pret ea et 3 : ee neies ae 
7 7 : een 
err 


oe ennnee 
TE te a 
eee EES 


Cente 
ia abel 


ap . set 
eee ne eer 0 
Be ee roa 


4m np ” 


